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1. INTRODUCTION 


The purpose of program VIEW is to compute view factors between specified 
surfaces and to be compatible with level 15. 5 of the NASTRAN 1 structural anal- 
ysis program. Program VIEW is a modification of a (finite element) view 
factor computation program called RAVFAC^. VIEW is designed to run on an 
IBM System/360 operating under OS (operating system), with a minimum region 
size of 110 K bytes. The actual computation of view factors is still performed 
exactly as it was in the original version of RAVFAC. In developing VIEW, 
RAVFAC was modified to satisfy the following compatibility requirements: 

• Accept finite element input which can also be used as input to 
NASTRAN. 

• Produce output (view factors) in a format which can be used as 
input to NASTRAN. 

• Follow NASTRAN program design so that in the future VIEW can be 
incorporated into NASTRAN as a subroutine. 

The VIEW program permits computation of the view factors between surfaces, 
taking into account the presence of any intermediate surfaces. VIEW also 
computes these view factors either by contour integration or by finite difference 
(double summation) methods. The first method is more accurate, but the 
second method is faster. Either method may be selected or a criterion may 
be specified which causes the program to select the best method based upon 
the geometry of the problem. 

Important features of the program include: 


• A restart capability, which protects a user against having to rerun 
an entire problem should a computer failure occur. 


NASTRAN is a structural analysis program developed by NASA at GSFC and 
is available from COSMIC. Level 15. 5 of NASTRAN includes thermal 
analysis capabilities. 

^RAVFAC was developed by J. K. Lovin and A. W. Lubkowitz, Lockheed 
Missile and Space Co,, Huntsville, Alabama, November 1969, Contract 
NAS8-30154. RAVFAC is available from COSMIC under No. MFS-21075. 
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• The ability to dynamically allocate available core space, thus 
allowing the user to request the amount of space in the computer 
required for his problem by using the region parameter on the job 
card. There is no maximum number of elements to which the user 
is limited, except that the computer's capacity may not be exceeded. 

• The ability to accept one or a combination of two input formats. A 
surface may be modeled by describing elements to program VIEW 
exactly as elements are described to RAVFAC or NASTRAN. 
Therefore, the VIEW program has RAVFAC-type and NASTRAN- 
type inputs. 

• The ability to run several problems in sequence in one job submis- 
sion. Each problem run is referred to as a "case". 

Programming information may be found in the "Programmer’s Manual for 
VIEW a Modification of the RAVFAC View Factor Program for Use with the 
NASTRAN Thermal Analyzer on IBM-360 Series Computers, ” GSFC Document 
No. X-322- 73-120, dated March 1973. 



2. VIEW PRIMER 


This section contains a general discussion on how to compute view factors, de- 
scribe the shape of elements, and describe the dimensions and locations of 
elements using RAVFAC- and NASTRAN-type data. The user should read this 
section to get an understanding of what type of information is to be input. This 
section does not explain how cards should be punched. The ordering of input 
and job control cards is also discussed. General information concerning the 
use of the VIEW program is presented in paragraph 2.6. 

2.1 COMPUTING VIEW FACTORS 

The computation of view factors is accomplished as follows: 

« The user constructs a finite element model of the surfaces for 
which he wishes to have the view factors computed. 

• The input data describes the location and shape of the finite 
elements relative to a basic coordinate system. 

• The program outputs the view factor from each element to each 
other element, taking into consideration the case when one or more 
elements block the view between two elements. 

An example of view factor computation is shown in Figure 2. 1. The user de- 
scribes the shape and location of elements 1, 2, and 3. The program computes 
view factors from El to E3, E3 to El, E3 to E2, E2 to E3, El to E2, and E2 
to El, taking into account that E3 partially shades El from E2. The equation 
used for computing these view factors may be found in Appendix E of the 
Programmer's Manual. (See page 1-2.) 



Figure 2.1. Example of View Factor Computation 
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2.2 AVAILABLE FINITE ELEMENT SHAPES 

The VIEW program allows the user to employ either NASTRAN-type input, 
RAVFAC-type, or a combination of both. 


Using only the NASTRAN-type input, the input data cards and punch output 
are compatible with the input to the NASTRAN program. Thus, a thermal 
analysis may be made using both the NASTRAN and VIEW programs, with data 
cards common to both. The availability of RAVFAC-type input is advantageous 
in some cases because it allows for more element shapes than are available 
with NASTRAN-type data. To be completely flexible, VIEW allows the use of 
both types of data in making a single model. The element shapes which may 
be described to the program according to the type of data that is used are 
given in Table 2.1. 


Table 2. 1 
Element Shapes 


NASTRAN-Type Data 

RAVFAC-Type Data 

Rectangular plate 

Rectangular plate 

Circular plate 

Circular plate or portions 

Triangular plate 

Trapezoidal plate 

Cylindrical shell 

Cylindrical shell or portions 

Conical shell 

Conical shell or portions 

Quadrilateral plate 

Spherical shell or portions 


Circular parabolic shell or 
portions 


These elements are described and illustrated in detail in Sections 4 and 5. 

The user specifies an element shape on a NASTRAN-type data card called 
CHBDY (see paragraph 6.3). If RAVFAC-type data is employed, then the 
element shape is specified on a RAVFAC-type data card number 2 (see para- 
graph 2.4). 
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2.3 DESCRIBING ELEMENTS USING NASTRAN-TYPE DATA 


2.3.1 Input Format 

The NASTRAN-type data cards used by VIEW are a subset of the cards known 
as "bulk data" cards in the NASTRAN program documentation. The NASTRAN- 
type input to the VIEW program must correspond to the following ruless 

• Cards are divided into ten eight-column fields as shown in Figure 2.2. 



Figure 2,2. Fields of a NASTRAN-Type Data Card 


• A mnemonic must be punched in field 1 beginning in column 1 to 
identify the card type. The mnemonics are listed in Section 6 and 
are comprised of names such as GRID, CHBDY, and CORDIS. 

• Fields 2 through 9 are for data items. All data items must lie com- 
pletely within a designated field, have no imbedded blanks, and must 
be of the proper type, that is, blank, integer, real, or BCD, All real 
numbers, including zero, must contain a decimal point. A blank will 
be interpreted as a real zero or integer zero as appropriate. Real 
numbers may be encoded in various ways. For example, the real 
number 7. 0 may be encoded as 7. 0, . 7E1, 0. 7 + 1, 70. -1, .70 + 1. 
Only the mnemonic in field 1 need be left- justified. 

• Normally field 10 is reserved for optional user identification. How- 
ever, the data required for a few types of cards do not always fit in 
fields 2 through 9. In this case, the remaining data is punched on a 
continuation card. In continuation cards, field 10, except column 73 
which is not referenced, is used in conjunction with field 1 of the 
continuation card as an identifier and hence must contain a unique 
entry. The continuation card must contain the symbol + in column 1, 
followed by the same seven characters that appeared in columns 

74 through 80 of field 10 of the card that is being continued. The con- 
tinuation card need not physically follow the card it continues 
in the NASTRAN-type data deck; however, run time will be reduced if 
it does. Figure 2.3 shows an example of a CHBDY card and its 
continuation card. 
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12 3456789 10 


Card 1 


Continuation 
of card 1 


Figure 2.3. Continuation of a NASTRAN CHBDY Card 

• Any NASTRAN-type data card which is mispunched in field 1 is ignored 
by VIEW. A message is printed out warning the user that "x" input 
cards were not recognized. The user should always check this 
message. 

2.3.2 Dimensions and Locations of Elements 

There are only five basic element shapes which may be described by using 
NASTRAN-type data. Each of these elements has a NASTRAN name as shown 
in Table 2.2. 


s 

CHBDY 

TTl 

X 

X 

X 

X 

X 

X 

X 

+ABC 

S' 

+ABC 

X 

X 

X 

X 

X 

X 

X 

X 



Table 2.2 

NASTRAN Element Names 


Element Shape 

NASTRAN Name 

Circular plate 

POINT element 

Rectangular plate 

LINE element 

Conical or cylindrical shell 

REV element 

Triangular plate 

ARE A3 element 

Quadrilateral plate 

AREA4 element 


Describing the dimension and location of these five elements is done by using 
grid points. Grid points are points in the three-dimensional space whose 
coordinates are specified on a NASTRAN-type data card called GRID (see 
paragraph 6. 3). The coordinates of a grid point may be specified relative to 
the basic rectangular coordinate system, or they may be specified relative to 
a local coordinate system which in turn is related to the basic coordinate 
system by use of NASTRAN-type data CORD cards. The local coordinate 
system may be either rectangular, spherical, or cylindrical (see paragraph 
6.3). 
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Example— The following grid points are located in space relative to a 
basic rectangular coordinate system and described to the program by the 
data cards shown in Figure 2.4. 



s 

GRID 

7 


1. 

1. 

1. 





GRID 

33 


i. 

1. 

0. 






Figure 2.4 Basic Coordinate System and Two Grid Cards 

The name and identification number of an element is punched on a NASTEAN- 
type data card called CHBDY. On this card the user also punches the grid 
point identification numbers of the grid points which determine the dimension 
and location of element 12. The specific format of the CHBDY card can be 
found in paragraph 6.3. 

Example— The user defines the size and location of a quadrilateral plate 
element he wishes to designate as number 12 by using the CHBDY card 
shown in Figure 2.5. 

The location of the grid points are specified on GRID cards. The grid points 
listed on the CHBDY card must be planar and must be listed sequentially. 

(The direction of ordering determines the active side of the element, as will 
be explained later . ) 

Example— The same comments in the previous example apply for the 
triangular element as shown in Figure 2. 6 but only three grid points need 
be on the CHBDY card. 
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The remaining elements cannot be completely defined using only grid points. 
Other information must appear on the CHBDY card. Paragraph 6.3 explains 
in detail how to specify all the elements using NASTRAN-type data. 

2.4 DESCRIBING ELEMENTS USING RAVFAC-TYPE DATA 

2.4.1 Input Format 

RAVFAC-type data allows the user to specify a surface shape which is divided 
into a mesh of elements by the program. (View factors are output on an 
element- to- element basis.) Thus, by specifying a surface shape, its location, 
dimensions, and an m-by-n mesh size, the user actually describes mn elements 
simultaneously. This saves the trouble of specifying mn separate shapes, 
locations, and dimensions. 

Three RAVFAC-type data cards are required to define each surface: cards 
number 1, 2, and 3. In card number 1 the user specifies the identification 
number of the surface, any desired comment, plus shading information. 

(Shading information is explained in Section 3. 1). Card 2 is used to define 
the surface shape, element mesh size and surface dimensions relative to a 
surface coordinate system. Card 3 is used to relate the surface coordinate 
system to the basic coordinate system. The specific format of each card is 
given in Section 6.4. 

Example— Figure 2.7 shows a RAVFAC-type data rectangular surface for 
which a two-by-three element mesh was specified. The program recognizes 
six adjoining rectangular elements, numbers them as described in the last 
part of Section 5, and then calculates the view factors to and from each one. 
A surface may be made one element by specifying a one-by-one mesh. 







RECTANGULAR 

ELEMENT 


RECTANGULAR SURFACE 


Figure 2.7. RAV FAC-Type Rectangular Surface Divided into Six Elements 
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2.4.2 Dimension and Location of Surfaces 


As discussed in paragraph 2.4. 1, the user specifies surfaces and element 
meshes on each surface. Knowing how an element mesh is defined over a 
surface, the user knows where each element is located once he has specified 
the location and dimensions of a surface. (See Section 5.) This paragraph will 
be concerned only with the question of how surfaces are specified. 


Each of the seven different RAVFAC surfaces has a surface coordinate system 
associated with it. Drawings of this concept can be seen in Section 5. Each 
surface is located in space by relating the surface coordinate system to the 
basic coordinate system of the model, or by relating it to a local coordinate 
system which in turn must be related to the basic coordinate system. 



Figure 2.8. Relationship of Basic and Surface Coordinate System 

Figure 2.8 shows the distances RX, RY, and RZ between the origins of the two 
coordinate systems. These distances, as well as PHI, PSI, and OMEGA, the 
angles which represent the rotation of the basic (or local) coordinate system 
into the surface coordinate system, must be specified on RAVFAC-type data 
card 3. The rotations must be listed in the order: 
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PHI = rotation of Y toward X about Z in degrees 
PSI = rotation of X toward Z about Y in degrees 
OMEGA = rotation of Y toward Z about X in degrees 

The size of a surface is specified by using the five parameters ALPHA, BMIN, 
BMAX, GMIN, and GMAX. These parameters define different surfaces in 
different ways as Figures 2. 9 and 2. 10 show. 


z 



z 



Figure 2. 10. Specification of a Cylindrical Surface 
See Section 5 for specification details of all surfaces. 
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2. 5 ORDERING INPUT AND JOB CONTROL CARDS 

Only a brief discussion on this subject will be made here. Detailed explanations 
are provided in paragraphs 6.1, 6.2, and 6.5. 

The set of input cards, of either one or both data types, which allows the pro- 
gram to compute the view factors of a given model is called a ’’case. " The. 
program allows the user to input several cases back-to-back, so that once 
the view factors for one case are found the program proceeds to the next case 
(Figure 2.11). 



Figure 2. 11. Data Deck Organization by Cases 

The cards which comprise the Job Control Cards deck in Figure 2. 11 are 
described in paragraph 6. 5. Each case is comprised of either five or seven 
case control cards and the NASTRAN- and/or RAVF AC-type input data. The 
format of the control cards is explained in paragraph 6.2. 

If a case has both types of input data, it must be ordered as shown in Figure 
2. 12. If a case has only one type of input data, there would be only that type 
of data and its control cards. A detailed explanation is given in paragraph 6. 1. 

2.6 GENERAL 

Although the method for describing the size and location of elements using 
either type of input data has been discussed, the program requires some 
additional information. All element shapes are two-dimensional. Although 
each element has two sides, the program will only compute view factors to and 
from one side. Therefore, the user must specify which side is active, that is. 
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Figure 2. 12. Deck Organization of More Than One Case 
Utilizing Both Types of Input Data 

to and from which side the view factors are to be computed. (The nonactive 
side can still cause shading. ) 

If the user is employing RAVFAC-type data, specify the active side on card 
type number 2 when indicating the surface shape using the ILK parameter. 
Refer to paragraph 6.4.2. For planar surfaces, a plus indicates the + Z 
side is active and a minus indicates the -Z side is active (Figure 2. 13). 
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For nonplanar surfaces, a plus indicates the outside is active while a minus 
indicates the inside is active (Figure 2. 14). 



Figure 2. 14. Indicating the Active Side of a RAVFAC Parabolic Element 


If NASTRAN-type data is being employed, indicate which side of an element is 
active on a CHBDY card. For the REV, AREA3, and AREA4 elements, the 
active side is specified by the ordering of the grid points. (Refer to the 
description of CHBDY cards in paragraph 6.3. ) For POINT and LINE elements, 
the active side is specified by the V vector on the CHBDY card (Figure 2.15). 



Figure 2. 15. Direction of V Indicates Active Side 


In Figure 2. 16 if grid point numbers are listed on the CHBDY and in the order 
1-2-3, 2-3-1, or 3-1-2, the active side is the one shown. A counter-clockwise 
ordering makes the other side active. Thus, for ARE A3 and AREA4 elements 
the righthand rule gives the active side .. 
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Figure 2. 16. Order of Grid Points Indicate Active Side 

A time-saving feature of the program should be emphasized. The program is 
capable of detecting whether one element causes shading between two other 
elements, as shown in Figure 2. 1. However, a significant amount of computer 
time can be saved if the user tells the program which surfaces can never cause 
shading and/or never be shaded. 

Using RAV FAC- type data, this specification is made on a card type number 1 
(see paragraph 6.4 for details). Using NASTRAN-type input data, this specifi- 
cation is made on a $VIEW card. The number of the $VIEW card is given in 
field 9 of the CHBDY card (see paragraph 6.3 for further details). 

The program also has the facility for creating an n-by-m mesh of subelements 
on a single element. View factors are calculated between subelements and 
internally recombined to print out the view factors between elements. Al- 
though this lengthens the run time it significantly increases the accuracy of 
answers, especially when shading occurs between two elements. This sub- 
element mesh size is specified on a $VIEW card if NASTRAN-type data is 
being used. It is specified on card type number 2 if RAVFAC-type data is 
used. At least a two-by-two subelement mesh is recommended on all elements 
until the user becomes quite familiar with the program (see Section 3 and 
paragraphs 6.3 and 6.4 for further details). A subelement mesh covers an 
element in the same manner that an element mesh covers a surface. It is 
recommended that the user look over the problems given in Appendixes A and 
B of this document. 

VIEW also has a restart capability. This means that a continuous record can 
be made on a magnetic tape while view factor computations are being made. 
Should some machine disaster cause termination before all computations 
have been completed, then the user may employ the magnetic tape to pick up 
his job from the stopping point. This saves the time and cost of having to 
restart a job from the beginning. (See paragraphs 6. 2 and 6. 5 for further 
details.) 
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The remainder of this document contains detailed information about each of 
the topics covered in this "primer” section. 
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3. MODELING 


The term modeling, as in modeling of a structure, is defined to mean the de- 
scription of a given structure to program VIEW, by approximating the structure 
using any combination of the following seven element shapes, and by specifying 
their size and location in space. 

• Triangular plate. 

• Quadrilateral plate. 

• Circular plate. 

• Cylindrical shell . 

• Conical shell. 

• Spherical shell. 

• Circular parabolic shell. 


NOTE 

The last two elements can be described only by using 
RAV FAC -type formatted input. 

As discussed in Section 2, element shapes may be described to the program by 
using either RAV FAC- or NASTRAN-type input. It is important to realize that 
program VIEW needs to know only element shapes, their geometrical location, 
and which side of the element is to be considered for view factor calculation. 
Hence it is possible to compute view factors between elements which do not form 
a valid structural or thermal model. 

3. 1 MODELING CONSIDERATIONS 

The user must specify which of the two sides of an element is to be considered 
active. That is, to and from which side are view factors to be computed. The 
inactive side will cause shading as well as the active side. 

Besides the possibility of modeling incorrectly in a structural or thermal sense, 
there are also some models which lead to erroneous view factor calculations 
for geometric reasons. This can best be explained by a brief description of 
how the program determines when a view factor between two elements is to be 
computed. 

In Figure 3. 1, V A is a vector normal to the centroid of element A, and points in 
the direction which the user specifies that element A can "see” and "be seen" 
by other elements; similarly for V g . A vector is then drawn from the centroid 
of B to the centroid of A, forming the angles a and 8 . Two conditions must be 
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satisfied simultaneously if there is to be a nonzero view factor between A and 
B: 

• -90° < 0 < +90° =*• cos 0 > 0 

• +90°< a < +180° or -90° > a > -180° =* cos a < 0 



Figure 3. 1. View Factor Computation Between Two Elements 


This is equivalent to stating that the active sides must see each other. Hence, 
whenever possible, the user should avoid having two elements meet anywhere 
except at a common edge. For example, if a structure included a wall meeting 
another wall to form a T-shape, it should not be modeled using three quadri- 
lateral elements, two (back-to-back) meeting the other at its center. Rather, 
four quadrilateral elements should be used, all meeting at a common edge. 


On the left in Figure 3. 2, element C, being right of center of element A, will 
have a view factor calculated from all of C to part of A. No view factor would 
be calculated between D and A. 



Figure 3. 2. Computing a Nonzero View Factor 
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In some cases it is impossible to satisfy this rule. For example, suppose the 
structure involves a plate whieh separates two cylinders, as shown in Figure 3. 3. 
The cylinder can and should be modeled as two separate cylinders which meet 
at the plate so that the previous rule is at least partially satisfied. The rule 
cannot be fully satisfied using the element shapes previously described. 



Figure 3. 3. Computing the View Factor for a Cylinder 

For such cases, the user should employ a feature of program VIEW which allows 
definition of a mesh of subelements * of any size on a particular element (Figure 
3.4). The view factors are calculated from these subelements and automatically 
recombined to produce a view factor for the element which they comprise. In 
fact, the program performs the "can I see” check described in Figure 3. 1 be- 
tween every set of two subelements . 



Figure 3.4 Subelement Mesh on a Plate 

By using this approach, inaccurate view factors may occur only between the 
cylinders and the subelements they touch. The remaining view factors between 
other subelements and the cylinders are found accurately. Thus, the smaller 


*For the user familiar with the RAVFAC program, the concept of subelements 
corresponds to elements in the RAVFAC documentation, while elements in this 
document correspond to RAVFAC nodes. In general, the terminology used in 
this document will follow that defined in NASTRAN, 
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the mesh the more accurate the view factors between the cylinders and elements, 
but the greater the run time. 

A subelement is shaded from another subelement if the line connecting the center 
of the two subelements is intersected by any part of at least one other element 
(Figure 3. 5). VIEW will automatically determine such situations if desired. 
However, the user may substantially reduce the problem run time if the elements 
which can shade and/or be shaded by other elements are specified. 



Figure 3. 5. Spherical Element B Shades Subelements A and C 

It is extremely important to realize that view factors are calculated on a sub- 
element-to- subelement basis and then recombined to print out a view factor from 
element to element. In Figure 3. 5, if the user specifies a one-by-one mesh of 
subelements on the elements A and C, a zero view factor would be calculated 
between A and C, although from the drawing this is not the true view factor. If 
the user specifies a finer mesh of subelements on elements A and C, as shown 
in Figure 3.6, then sphere B, which shades mostly subelement A2 from C2, 
would cause a zero view factor to be computed only between those two subele- 
ments. A nonzero view factor would be computed from each subelement of A to 
each subelement of C (except between A 2 and C2), and then recombined to give 
a view factor between elements A and C. The user should note from this exam- 
ple that again, the finer the mesh of subelements the more accurate the view 
factors. 



Figure 3. 6. Computing View Factor for a Fine Mesh of Subelements 




A subelement mesh must be specified for all nonplanar elements. The subele- 
ments for nonplanar elements are recognized by the program as flat rectangular 
plates and view factors are calculated from and to these plates. A cylinder 
modeled with an eight -by-four mesh, as shown on the left in Figure 3.7, is 
recognized by the program as shown on the right in Figure 3. 7. 



Figure 3.7. Cylinder Modeled With Eight-by-Four Mesh 

While a fine mesh of subelements will produce more accurate answers than a 
coarse mesh, it will also increase computer run time. When not needed, a fine 
mesh of subelements should be avoided. An example of such a case is Figure 3. 6 
if no element shaded A from C. Here a one-by-one mesh of subelements 
on elements A and C should be adequate. 

3. 2 COORDINATE SYSTEMS 

VIEW allows the user the convenience of specifying the location of different sur- 
faces in different coordinate systems. In describing elements using NASTRAN 
input, coordinate systems are described and used just as in the NASTRAN 
program. Similarly, in describing surfaces using RAVFAC input, coordinate 
systems are described and used just as in the RAVFAC program. 

Both types of input allow the user to specify local coordinate systems for con- 
venience. That is, the user is allowed to locate other coordinate systems in 
space by specifying their orientation and displacement from the basic coordinate 
system. By using NASTRAN data, the user has the added choice of selecting 
the type of local coordinate system he desires. It may be either rectangular, 
spherical, or cylindrical. 

A limitation imposed by VIEW which is not present in NASTRAN is that each 
local coordinate system must reference the basic coordinate system. One 
local coordinate system may not be specified relative to another local coordinate 
system. Also, the program will allow the use of, at most, forty NASTRAN 
local coordinate systems of type CORD 1 and forty of type CORD 2. 
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Figure 3. 8 shows all the basic methods of describing the location of elements 
of a model. Each model must have one and only one basic coordinate system. 
For NASTRAN-type input each element of the model may be located in space by 
specifying the coordinates of bounding grid points relative to the basic coordinate 
system. For RAVFAC-type input, the user relates the location and orientation 
of an element coordinate system to the basic coordinate system. 


Z 



Figure 3.8. Various Methods of Describing Locations of Elements of a Model 

As an extra feature of the program, the user may relate grid points or surface 
coordinate systems (depending on which type of input he is using) to one or more 
local coordinate systems. He then specifies the location and orientation of the 
local coordinate system to the basic coordinate system. 
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4. NA3TRAN ELEMENTS 


There are five element shapes which may be described to VIEW by using NAS- 
TRAN data (Table 4. 1). More shapes are available using RAVFAC data (see 
Section 5). Either side of an element, but not both sides, may be designated as 
active, that is, the side to and from which view factors are computed. Specifi- 
cation of the active side is described in paragraph 6. 3. 


Table 4. 1 

NASTRAN Elements 


Element Shape 

NASTRAN Name 

Circular plate 

POINT 

Rectangular plate 

LINE 

Cone or cylinder 

REV 

Triangular plate 

A REA 3 

Quadrilateral plate 

A REA 4 


There is one additional element shape which may be described to NASTRAN for 
heat transfer analysis purposes, but which is not available in VIEW, that is, the 
elliptic cylinder. A circular cylindrical shape is available by using the REV 
element. Figures 4. 1 through 4. 6 illustrate the element shapes (broad lines) 
and an example of their division into subelements (fine lines). 



Figure 4. 1. Circular Plate With Two-by-Four Mesh of Subelements 

In Figure 4. 1, NB is the division in circular sections while NG divides the disk 
into pie-shaped sections. 
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• LINE Element 



NB = 4 
NG = 3 


Figure 4. 2. Rectangular Plate With Four-by-Three Mesh of Subelements 

In Figure 4.2, NB is the division perpendicular to the line connecting the two 
grid points needed to give the location of the LINE element. NG is the division 
along the line. 

• REV Element 



Figure 4. 3. Cylinder With Two-by-Six Mesh of Subelements 

In Figure 4. 3, NB is the division along the axis of the element, while NG is the 
division around the surface. NASTRAN requires that grid points G1 and G2 be 
in the X-Z plane of the basic coordinate system. 
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e AREA 3 Element 


G3 



Figure 4.4. Triangular Plate With Seven-by-Three Mesh of Subelements 

In Figure 4.4, the dashed line is a perpendicular from G3 to the line connecting 
G1 and G2. NB is the division along this line. NG is the division along the 
line connecting G1 and G2. 


• AREA 4 Element 



Figure 4. 5. Nonrectangular Quadrilateral Plate With Specified Four-by-Three 

Mesh of Subelements 

The complexity of Figure 4. 5 is due to the fact that the element is divided into 
two triangular elements by connecting grid points G1 and G3. The subelement 
mesh size specified by the user (in this case NB-^ - 4 and NG^ = 3) is applied to 
triangle Gl, G2, G3, just an described for the AREA3 element. NG divisions 
are made along the line connecting Gl and G2, while NB divisions are made along 
a perpendicular to this line. However, for triangle G3, G4, Gl, the subelement 
mesh size is chosen proportionate to that specified for the first triangle in the 
following manner: 


NGjl ng 2 

length (Gl, G2) = length (G3, G4) 

NB X NB 2 

length (G3, X) ” length (Gl, Y) 
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If either NB 2 or NGg are calculated to be zero, then both NBg ancj NG2 are set to 
one and a warning message is printed. Note that in counting the number of 
elements in the problem, each nonrectangular AREA4 element should be 
counted as two elements. If the AREA4 elements is rectangular, then it is not 
divided into two triangular elements by the program. 



NB = 3 
NG = 4 


Figure 4. 6. Mesh Formed if AREA4 Element is Rectangular 


In Figure 4. 6, NB divisions are again made along the line connecting G1 and G2, 
while NG divisions are made along the line connecting G2 and G3. In this case, 
the AREA4 element is counted as one element, as it is not broken into two 
triangular elements. 
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5. RAVFAC SURFACES 


The VIEW program also accepts input in the same format as the RAVFAC pro- 
gram. (See Section 2. ) The input format is described in detail in paragraph 
6.4. While NASTRAN-type input allows the user to describe elements and the 
mesh of subelements into which each element is divided, RAVFAC-type input 
is somewhat more general. 

The user may describe a surface comprised of a mesh of elements, each of 
which is comprised of a mesh of subelements. Hence, where with NASTRAN 
input the user might model a wall with 100 rectangular elements (10 along the 
base and 10 along an edge), with RAVFAC-type input the user may form the 
same model by specifying one surface comprised of a ten-by-ten mesh of 
elements. View factors are still output on an element- to-element basis. 

Another advantage of RAVFAC-type input is that more surface shapes may be 
described than by using NASTRAN-type input. These surface types are: 

• Rectangular plates. 

• Circular plate or section thereof. 

• Trapezoidal plate. 

• Cylindrical shell or portion thereof. 

• Conical shell or portion thereof. 

• Spherical shell or portion thereof. 

• Circular parabolic shell or portion thereof. 

Each element is described by specifying its location and active side relative 
to the element coordinate system as shown in Figures 5. 1 through 5.8. The 
element coordinate system is then related to the local or basic coordinate 
system. Refer to Section 3 and paragraph 6.4 for further details on specifi- 
cation of the active side. 
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Figure 5.4. Surface Type ±4 Cylinder 



Figure 5. 5 Surface Type ±5 Cone 
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Each surface may be divided into NVB elements in the /? direction and NVG 
elements in the y direction. In turn each element may be divided into NB 
subelements in the (3 direction and NG subelements in the y direction. 



Figure 5.8. Four- By-Three Element Mesh with Three-By- 
Four Subelement Mesh 

On element n + 5 of Figure 5. 8 is an example of a three-by-four subelement 
mesh. The user may label the surface any number n, but the program will 
number the elements as indicated. 
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6, INPUT 


6. 1 ORGANIZATION 

Several view problems may be run at one time. Each problem is referred to 
as a case. A typical input data deck would appear as shown in Figure 6. 1. 



— — CASE 2 DATA 
CASE 1 DATA 


Figure 6.1. Input Data Deck 

Any case must always begin with a title card, after which must come a case 
control card followed by a format- type card declaring NASTRAN input or 
RAVFAC input is to follow. If NASTRAN input is used, after the NASTRAN 
data must come an ENDDATA card (Figure 6.2). If RAVFAC data is used, 
the end of data is signified by inserting a card with -1 in columns 4 and 5 
(Figure 6.3). If the user desires to run a problem using both formats, 
NASTRAN data must appear before the RAVFAC data (Figure 6.4), Every 
case must end with an ENDCASE card so that the program may distinguish 
between cases. 

Figures 6.2 through 6.4 describe the only ways that the control cards com- 
prising a case may be organized. Any data card may be punched on either 
the 026 or 029 key punch. 

6.2 CONTROL CARDS 

The control cards portion of a case is described in the following subparagraphs. 
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ENDCASE CARD 
ENDDATA CARD 


NASTRAN INPUT 
FORMAT TYPE CARD 

CASE CONTROL CARD 
TITLE CARD 


Figure 6. 2. Organization of a Case Consisting Only of NASTRAN Data 



ENDCASE CARD 
SIGNALS END OF RAVFAC INPUT 


RAVFAC INPUT 


FORMAT TYPE CARD 


CASE CONTROL CARD 


TITLE CARD 


Figure 6. 3. Organization of a Case Consisting Only of RAVFAC Data 
6.2.1 Title Card 

This is the first card in a case. It must have TITLE punched in columns 1 to 5. 
Beginning in column 9, the user may punch any 71 characters, and those 
characters are printed at the beginning of the output to identify the case. See 
sample output in appendices. 
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RAVFAC 
INPUT 


NASTRAN 
INPUT 


f RAVFAC INP UT 



ENDCASE CARD 
SIGNALS END OF RAVFAC INPUT 


FORMAT TYPE CARD 


ENDDATA CARD 


FORMAT TYPE CARD 
CASE CONTROL CARD 
TITLE CARD 


Figure 6.4. Organization of a Case Consisting of Both NASTRAN 

and RAVFAC Data 


6.2.2 Case Control Card 

This is the card which immediately follows the title card and is used to pass 
information to the program regarding restart, type of integration to be used, 
and various other features as described in the following format description. 
The format of this card is as follows: 

Columns 1 to 5 — User must punch CASE=. This identifies the 

card as a case control card. 

Columns 9 to 16 — Enter integer anywhere in columns 9 to 16, which is 

the same number of minutes CPU time shown on the 
JCL job card. Punch H if the time estimate on the 
job card is one-half minute. Defaults to 100,000, 
but will not affect results unless job times out, in 
which case all available output may not be printed 
out. If the proper CPU time is punched in this 
field and the job times out then all computed view 
factors will be printed out. 

Columns 17 to 24 — An integer less than or equal to zero in these 

columns indicates that no restart tape is to be 
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produced. Any integer greater than zero indicates 
that a restart tape will be produced on unit 2. In 
this case the user must supply the computer opera- 
tors with a tape and specify the number of the tape 
on the EXEC card in the JCL (paragraph 6.5). 

Columns 25 to 32 — If any integer greater than zero appears in columns, 
then all view factor computations are made using the 
finite difference technique. An integer less than 
zero causes contour integration to be used in com- 
puting all view factors. If zero is punched in columns 
25 through 32, or if they are left blank, the program 
selects which of the two methods it will use between 
two elements, based on the parameter RMAX speci- 
fied in columns 49 through 56. Consider the situa- 
tion in Figure 6.5. 



d 


rs 



Figure 6, 5. Choosing a View Factor Computation Method 

The finite difference method becomes inaccurate 

o 

when d rs >A g /0. 1 where A g is the area of a 
subelement s, and d rs is the distance between 
two subelements for which view factors are being 
computed. The solution to this problem is not to 
use contour integration for all calculations be- 
cause it requires more computer time. By speci- 
fying a zero in any one of columns 25 through 32, 
the program selects the most efficient techniques 
between two elements by checking whether: 

A 

!_ > RMAX 


where RMAX is specified in columns 49 through 56, 
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If this relation is satisfied, contour integration is 
used. If not, finite difference is used. This check 
requires more time than if the user had requested 
the finite difference technique, but less time than 
if contour integration was requested. 

Columns 33 to 40 — Any integer greater than zero causes the subelement 
mesh size for each element to be set to one-by-one, 
overriding information given on the input data. For 
large problems a short checkout run of the data 
validity should be made. Using this option to set 
a one-by-one subelement mesh for each element will 
allow for the problem to be processed in a relatively 
short time. A zero, blank, or negative integer 
tells the program that the mesh sizes will be speci- 
fied on the input data. This should be used once 
the user is sure of the model. 

Columns 41 to 48-— Any integer less than zero indicates that all shading 
is to be neglected. This will save computer time, 
but the user must be sure that no shading occurs 
in the model. An integer greater than or equal to 
zero tells the program to consider shading, using 
the information which describes each surface in 
the input data. There the user can specify whether 
or not a surface can shade any other surfaces or 
be shaded by any other surfaces. 

Columns 49 to 56 — This is where the user specifies the floating point 
number RMAX. If the columns are left blank, 

RMAX is set to 0. 1 by the program. A value of 
. 01 will provide results as accurate as usually 
required. A value greater than . 1 will probably 
produce view factors which are too large. If the 
user selects contour integration or finite differ- 
ence integration in columns 17 to 24 rather than 
allowing the program to choose the method, then 
any number in these columns is ignored. 
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6.2.3 Format Type Card 


This is the card which always appears after the case control card and some- 
times after an ENDDATA card (see Figure 6.4). Its purpose is to identify 
the format type of the data which follows. Punch BEGIN BULK in columns 1 
to 10 to indicate NASTRAN-type input will follow. Punch RAVFAC DATA in 
columns 1 to 11 to indicate that RAVFAC -type data will follow. 

6.2.4 ENDDATA Card 

This card immediately follows the last card input in NASTRAN format. It 
signifies the end of NASTRAN-formatted input for the case. The word 
ENDDATA is punched beginning in column 1. 

6. 2. 5 -1 Card 

This card immediately follows the last card input in RAVFAC format. It 
signifies the end of the RAV FAC-formatted input for the case. Punch -1 in 
columns 4 and 5. 

6.2.6 ENDCASE Card 

This card follows either an ENDDATA card or a -1 card. It signifies the end 
of the case. The word ENDCASE is punched beginning in column 1. 

6.2.7 -1 Card (For Local Coordinate Systems) 

This card follows the last RAVFAC local coordinate system card (paragraph 
6.4.4). It is used only if RAVFAC local coordinate system cards are used. 
Punch -1 in columns 9 and 10. 
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6.3 NASTRAN- FORMATTED INPUT DATA CARDS 


Once the program has read a format-type card declaring BEGIN BULK, it 
knows that it will be reading input data cards in a format which is also accepta- 
ble to the NASTRAN program. At this point, VIEW is similar to NASTRAN. 
Some of the read subroutines used in VIEW were taken directly from the 
NASTRAN program. 

The NASTRAN input data cards may be arranged in any order, but processing 
time will be decreased if continuation cards follow directly behind the cards 
which they continue. The free field format feature is also available in VIEW. 
That is, integer, decimal, alphanumeric, and exponential numbers may be 
punched starting in any of the columns of fields 2 through 10. * Alphanumeric 
characters must be left-justified in field 1. The large field option available in 
NASTRAN for the input of double precision is not available in VIEW. Cards 
may be punched on an 026 or 029 keypunch. See paragraph 2. 3. 1 for further 
details. 

Once VIEW realizes it is to read NASTRAN- formatted data, it will accept any 
type of card without causing an error until it reaches an ENDDATA card. 
However, VIEW will only use and print out the following types of cards: 

• CHBDY 

• PHBDY 

• GRID 

• CORD1R, S, C 

• CORD2R, S, C 
o GRDSET 

• $VIEW 


NOTE 

If any other types of cards are read a message is 
printed out on the last page of output stating: 

X LOGICAL CARDS DETECTED OF TYPES NOT 
RECOGNIZED BY VIEW. 


*A NASTRAN input data card is divided into ten fields, each consisting of eight 
columns. A continuation to a card is organized in the same fashion and its 
fields are numbered 11 to 20. 
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6.3.1 CHBDY Input Data Card 


This card defines an element, numbers it, references, and gives data concern- 
ing its shape and geometric location. Figure 6. 6 is the format and an example 
of a CHBDY card. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CHBDY 

EID 

PID 

TYPE 

Gl 

G2 

G3 

G4 

IVIEW 

+abc 

CHBDY 

413 

29 

LINE 

107 

108 



14 

+BD27 


12 3 4567 89 10 


+abc 

GA1 


GA3 

GA4 

VI 

V2 

V3 



+BD27 


1 



1.00 

0.0 





Figure 6. 6. Format and Example CHBDY Card 
The contents of each CHBDY input data card are as follows: 



Field 

Contents 

i. 

CHBDY 

CHBDY should be punched left- justified to identify 



the card type. 

2. 

EID 

Element identification number (Integer >0) 

3. 

PID 

Integer > 0 or blank. Only if the element being 


specified is either a POINT or LINE element 
must the user reference a PHBDY card. The 
PHBDY card gives the area of a POINT element 
or the width of a LINE element. 

4. TYPE Defines the shape of the element. Punch one of 

the names POINT, LINE, REV, ARE A3, AREA4. 

5 to 8. Integer > 0 or blank. Grid point identification 

Gl, G2, G3, G4 numbers. The boundaries of the element will be 

determined by the grid points listed here. 
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Field 


Contents 


9. IVIEW Integer >0. $VTEW card identification number. 

The element will have the characteristics given 
by the $VIEW card with this number. 

This card must be continued only if it defines a POINT or LINE element. If 
one of these two elements is defined and the card is not continued, the element 
will be ignored by the program and will not appear in the view factor calcula- 
tions. If the card is to be continued, the following fields must be filled in: 

Field Contents 


10, 11. +abc 


12 to 15. 

16 to 18. 

VI, V2, V3 


Field 10 of the CHBDY card (except column 73 
which is not referenced) is used in conjunction 
with field 11 of a continuation card as an identi- 
fier, and hence must contain a unique entry. The 
continuation card must contain the symbol + 
in column 1, followed by the same seven charac- 
ters that appeared in columns 74 to 80 of field 10 
of the card that is being continued. See para- 
graph 2.3. 1 for further details. 

Anything punched in these fields is ignored. 

In these fields, VI, V2, and V3 are the coordi- 
nates of the vector in the coordinate system 
referenced in field 7 of grid card Gl. The 
direction of the vector orients the elements and 
indicates the active side as described on the 
next page . 


The following are descriptions of the element types. As seen in Figure 6. 7, 
the POINT element shape is a flat disk with the center at grid point Gl. Its 
active side and orientation in space are given by the vector, specified in 
fields 16, 17, and 18 of the CHBDY continuation card. The ^ is given in the 
coordinate system referenced in field 7 of grid card Gl. Specify the disk area 
by using a PHBDY card. 



Figure 6.7. POINT Element — Flat Disk Shape 
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The LINE element shape, as seen in Figure 6.8, is a flat rectangular plate. 
Grid points G1 and G2 are located at the center of two opposite sides. The 
specified vector V and the vector T form a plane. The LINE element is 
oriented perpendicular to this plane. The active^side^of the element is given by 
the normalized vector n = (TX(VXT) )/|TX(VXT)| . Again, V is given in 
the coordinate system referenced in field 7 of grid card Gl. Specify the width 
perpendicular to T by using a PHBDY card. 



Figure 6.8. LINE Element — Flat Rectangular Plate 

As seen in Figure 6. 9, the REV element shape is a three-dimensional cone or 
cylinder. The grid points Gl and G2 must be in the X-Z plane of the basic 
coordinate system. A vector T from Gl to G2 connects them. The element is 
generated by revolving T about the Z axis 360 degrees. The active side is 
given by the normalized vector n = (!y X T)/ | ey X t| . There is no need to 
specify V or a surface area. 


z 



Figure 6.9. REV Element — Three-Dimensional Cone 
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The AREA3 element shape is a flat triangular plate, as shown in Figure 6. 10. 
It is oriented in space by the locations of grid points Gl, G2, and G3. Its 
active side is specified by the normalized vector n = (T 12 X Tj 3 )/ |t 1 2 X T-^ 3 )| 
(that is^the right-hand rule is applied to Gl, G2, G3). There is no need to 
specify V or a surface area. 


G3 



Figure 6.10. AREA3 Element — Fla.t Triangular Plate 


As seen in Figure 6. 11, the AREA4 element shape is a flat quadrilateral plate. 
It is oriented in space by the locations of Gl, G2, G3, and G4. The grid points 
must all lie in on e plane. They must be entered so that the lines G1G2, G2G3, 
G3G4, and G4G1 do not intersect except at grid points. The active side is 
determined^ in the same manner as for the AREA3 element. There is no need 
to specify V or a surface area. 



Figure 6.11. AREA4 Element — Flat Quadrilateral Plate 
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6.3.2 CORD 1C Input Data Card 


The CORD 1C card defines a cylindrical coordinate system by reference to 
three grid points. These points must be defined in the basic coordinate sys- 
tem. The first point is the origin, the second lies on the Z-axis, and the 
third lies in the plane of the azimuthal origin. Figure 6. 12 shows a cylindrical 
coordinate system and Figure 6. 13 is the format and an example of a CORD 1C 
card. 


z 



Figure 6. 12. Cylindrical Coordinate System 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

COR D1C 

C1D 

G1 

G2 

G3 

CID 

G1 

G2 

G3 


COR D1C 

3 

16 

32 

19 







Figure 6.13. Format and Example CORD 1C Card 
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The contents of each CORD 1C input data card are as follows: 


Field 


Contents 


1. CORD 1C Punch CORD 1C beginning in column 1 

2. CID Coordinate system identification number 

(Integer > 0) 

3 to 5. Gl, G2, G3 Grid point identification numbers (Integer > 0); 

G1 * G2 * G3) 

Coordinate system identification numbers on all CORD1C cards must be 
unique. The three points Gl, G2, G3 must be noncollinear and in the basic 
coordinate system. The location of a grid point (P in Figure 6. 12) in this 
coordinate system is given by R, 0, TJ where G is measured in degrees. 
One or two coordinate systems may be defined on a single card. 
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6.3.3 CORD1R Input Data Card 


The CORD 1R card defines a rectangular coordinate system by reference to 
three grid points. These points must be defined in the basic coordinate sys- 
tem. The first point is the origin, the second lies anywhere on the +Z-axis, 
and the third lies anywhere in the X-Z plane. Figure 6. 14 shows a rectan- 
gular coordinate system and Figure 6. 15 is the format and an example of a 
CORD1R card. 


z 



Figure 6.14. Rectangular Coordinate System 
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Figure 6. 15. Format and Example CORD1R Card 
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The contents of each CORD lit input data card are as follows: 


Field 


Contents 


1. CORD1R Punch CORD1R beginning in column 1 

2. CID Coordinate system identification number 

(Integer > 0) 

3 to 5. Gl, G2, G3 Grid point identification numbers (Integer >0); 

G1 * G2 G3) 

Coordinate system identification numbers on all CORD1R cards must be 
unique. The three points Gl, G2, G3 must be noncollinear and in the basic 
coordinate system. The location of a grid point (P in Figure 6. 14) in this 
coordinate system is given by X', Y' t Z' . One or two coordinate systems 
may be defined on a single CORD1R card. 
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6.3,4 CORDIS Input Data Card 


The CORDIS card defines a spherical coordinate system by reference to three 
grid points. These points must be defined in the basic coordinate systems. 

The first point is the origin, the second lies on the +z-axis, and the third lies 
in the plane of the azimuthal origin. Figure 6. 16 shows a spherical coordinate 
system and Figure 6. 17 is the format and an example of a CORDIS card. 
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Figure 6. 16. Spherical Coordinate System 
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Figure 6. 17. Format and Example CORDIS Card 
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The contents of each CORDIS input data card are as follows: 

Field Contents 

1. CORDIS Punch CORDIS beginning in column 1. 

2. CID Coordinate system indentification number 

(Integer > 0) 

3 to 5. Gl, G2, G3 Grid point identification numbers (Integer > 0; 

G1 5* G2 ^ G3) 

Coordinate system identification numbers on all CORDIS cards must be 
unique. The three points Gl, G2, G3 must be noncollinear and in the basic 
coordinate system. The location of a grid point (P in Figure 6. 16) in this 
coordinate system is given by R, 9, 4* where 0 and are measured in 
degrees. One or two coordinate systems may be defined on a single 
CORDIS card. 
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6.3.5 C0RD2C Input Data Card 


The C0RD2C card defines a cylindrical coordinate system by reference to the 
coordinates of three points. The first defines the origin, the second defines 
the direction of the +Z-axis, and the third lies in the plane of the azimuthal 
origin. Figure 6. 18 shows a cylindrical coordinate system and Figure 6. 19 
is the format and an example of a CORD2C card. 
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Figure 6.18. Cylindrical Coordinate System 
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Figure 6. 19. Format and Example of CORD2C Card 


CARO 1 


CARD 2 
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The contents of each CORD2C input data card are as follows: 


Field 


Contents 


1 . 

CORD2C 

2. 

CID 

3. 

RID 

4-6. 

Al, A2, A3 

7-9. 

Bl, B2, B3 

12-14. 

Cl, C2, C3 

10, 11. 

+ABC 


Punch CORD2C beginning in column 1 
Coordinate system identification number 
(Integer > 0) 

Blank 

Coordinates of three points in the basic 
coordinate system (Real) 

Unique continuation symbol. (See 2.3. 1) 


The continuation card must be present. The three points (Al, A2, A3), (Bl, 
B2, B3), (Cl, C2, C3) must be unique and noncollinear. Noncollinearity is 
checked by the geometry processor. The points must be in the basic coordi- 
nate system. Coordinate system identification numbers on all CORD2C cards 
must be unique. The location of a grid point (P in Figure 6. 18) in this 
coordinate system is given by R, 9, Z' where 0 is measured in degrees. 
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6.3.6 CORD2R Input Data Card 


The CORD2R card defines a rectangular coordinate system by reference to 
the coordinates of three points. The first point defines the origin, the second 
defines the direction of the +Z-axis, and the third defines a vector which, with 
the Z-axis, defines the X-Z plane. Figure 6.20 shows a rectangular coordi- 
nate system and Figure 6. 21 is the format and an example of a CORD2R card. 
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Figure 6.20. Rectangular Coordinate System 
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Figure 6.21. Format and Example CORD2R Card 


CARD 1 


CARD 2 
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The contents of each CORD2R input data card are as follows: 


Field 


Contents 


1. 

CORD2R 

2. 

CID 


3. 

RID 


4-6. 

Al, A2, 

A3 

7-9. 

Bl, B2, 

B3 

12-14. 

Cl, C2, 

C3 

10, 11. 

+ABC 



Punch CORD2R beginning in column 1 
Coordinate system identification number 
(Integer > 0) 

Blank 

Coordinates of three points in the basic 
coordinate system (Real) 

Unique continuation symbol. (See 2. 3. 1) 


The continuation card must be present. The three points (Al, A2, A3), (Bl, 
B2, B3), (Cl, C2, C3) must be unique and noncollinear; noncollinearity is 
checked by the geometry processor. The points must be in the basic coordi- 
nate system. Coordinate system identification numbers on all CORD2R cards 
must be unique. The location of a grid point (P in Figure 6. 20) in this coor- 
dinate system is given by X', Y', Z'. 
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6.3.7 CORD2S Input Data Card 


The CORD2S card defines a spherical coordinate system by reference to the 
coordinates of three points. The first point defines the origin, the second de- 
fines the direction of the +Z-axis, and the third lies in the plane of the azi- 
muthal origin. Figure 6.22 shows a spherical coordinate system and Figure 
6. 23 is the format and an example of a CORD2S card. 
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Figure 6.22. Spherical Coordinate System 
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Figure 6.23. Format and Example CORD2S Card 
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The contents of each CORD2S input data card are as follows: 


Field 


Contents 


1 . 

CORD2S 

2. 

CID 

3. 

RID 

4-6. 

Al, A2, A3 

7-9. 

Bl, B2, B3 

12-14. 

Cl, C2, C3 

10, 11. 

+ABC 


Punch CORD2S beginning in column 1 
Coordinate system identification number 
(Integer > 0) 

Blank 

Coordinates of three points in the basic 
coordinate system (Real) 

Unique continuation symbol. (See 2.3.1) 


The continuation card must be present. The three points (Al, A2, A3), (Bl, 
B2, B3), (Cl, C2, C3) must be unique and noncollinear; noncollinearity is 
checked by the geometry processor. Coordinate system identification numbers 
on all CORD2S cards must be unique. The location of a grid point (P in 
Figure 6.22) in this coordinate system is given by R, 6, 4> where 6 and <t> 
are measured in degrees. 
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6.3.8 GRID Input Data Card 


GRID cards define the location of a geometric grid point of the model. Figure 
6. 24 is the format and an example of a GRID card. 


1 
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GRID 

ID 

CP 

XI 

X2 
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CD 
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GRID 
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6 
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0.0 

-4.5 






Figure 6. 24, Format and Example GRID Card 
The contents of each GRID input data card are as follows! 

Field Contents 


1. GRID 

2. ID 

3. CP 

4, 5, 6. XI, X2,X3 


GRID should be punched left-justified to identify the 
card type. 

Grid point identification number (Integer >0). 

Identification number of the coordinate system in 
which the location of the grid point is defined 
(Integer >0). 

Location of the grid point in coordinate system CP 
(any real number). 


7. CD Identification number of a coordinate system in 

which the V vector on a CHBDY card is defined. 
Used only for Point and LINE elements (Integer >0). 

8 to 10. Any numbers punched in these fields are ignored 

by VIEW. 


The meaning of XI, X2, and X3 depend on the type of coordinate system 
defined by CORD card CP. The system may be rectangular, cylindrical, or 
spherical coordinates. 
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6.3.9 GRDSET Input Data Card 

This card defines default options for fields 3, 7, and 8 of all grid cards. 
Figure 6. 25 is the format and an example of a GRDSET card. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

GRDSET 


CP 


@1 


CP 

PS 



GRDSET 


16 




32 





Figure 6.25. Format and Example GRDSET Card 


The contents of each GRDSET input data card are as follows: 

Field Contents 

1. GRDSET GRDSET should be punched left- justified to identify the card 
type. 

3. CP Identification number of the coordinate system in which the 

location of the grid point is defined (Integer >0). 

7 . CD Identification number of a coordinate system in which the V 

vector on a CHBDY card is defined (Integer >0). 

8. PS Not used by VIEW. 

The contents of fields 3, 7 and 8 of this card are assumed for the corresponding 
fields of any grid card whose fields 3, 7 and 8 are blank. Only one GRDSET 
card may be used for each case. 
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6.3.10 PHBDY Input Data Card 


The PHBDY card defines properties of elements specified on CHBDY cards, 
and is only needed for POINT and LINE elements. Figure 6. 26 is the format 
and an example of a PHBDY card. 
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PHBDY 
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Figure 6. 26. Format and Example PHBDY Card 
The contents of each PHBDY input data card are as follows: 



Field 

Contents 

i. 

PHBDY 

PHBDY should be punched left -justified to identify the card 



type. 

2. 

PID 

Property identification number (Integer >0). 

4. 

AF 

Area if referenced by a CHBDY card defining a POINT element. 


Width if referenced by a CHBDY card defining a LINE element. 
(Real ^0.0), 


Note: Any data punched in fields 3 and 5 to 10 
are ignored by VIEW. 
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6.3.11 $VIEW Input Data Card 


This card is used only by the VIEW program and is not recognized by the 
NASTRAN program. It is used to give element characteristics which are not 
found on the other NASTRAN cards. Figure 6. 27 is the format and an example 
of a $VIEW card. 
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Figure 6.27. Format and Example $VIEW Card 
The contents of each $VIEW input data card are as follows: 


Field 

1. $VIEW 

2. IVIEW 

3. KSHD 

4. KBSHD 

5. NB 

6. NG 

7. DISLIN 


Contents 


$VIEW should be punched left- justified to identify the card 
type. 

VIEW card identification number (Integer >0). 


Can shade flag (0 = NO, 1 = YES). 

Can be shaded flag' (0 = NO, 1 = YES). 


Subelement mesh size in the beta direction (Integer >0). 
Subelement mesh size in the gamma direction (Integer >0). 


Displacement of a surface perpendicular to the active side of 
the surface. Used only for the LINE element (Real). 


The shading flags should be used carefully. When the user can identify a sur- 
face which cannot cause shading between any other surfaces, he can save com- 
puter time by identifying it. If in doubt as to whether or not it can shade, it 
should be labeled "can shade." Similarly, the user can save computer time by 
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identifying surfaces which cannot be shaded. If in doubt, they should be labeled 
"can be shaded.” If the user wants the program to do all the work, all surfaces 
must be labeled KSHD =1, KBS HD =1. 

The DISLIN field should be used only for LINE elements. It locates the surface 
DISLIN units perpendicular to the place where the active side of the surface is 
located by the CHBDY card (Figure 6. 28). 


DISLIN 


Figure 6. 28. Location of LINE Elements For DISLIN < 0 
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6.4 RAV FAC- FORMATTED INPUT DATA CARDS 


The description of the geometry of a surface using RAV FAC-formatted input 
is totally different from that of an element using NASTRAN-formatted data. 
However, the information is essentially identical. 


When using RAVFAC data, the user needs three distinct cards to describe 
each surface. These sets of three cards may be entered in any order; however, 
the ordering of the three cards within a set is required to follow a certain 
rule (refer to Figure 6.29). 


L 


z 




£ 


■1 


L 


RAVFAC DATA 


1 PX 


3 
21— J 


2|fj 


|. IDENTIFIES END OF RAVFAC 
INPUT (SEE PARAGRAPH 6.2) 


DESCRIBE SECOND INPUT SURFACE 


DESCRIBE FIRST INPUT SURFACE 

IDENTIFIES FORMAT OF DATA 
TO FOLLOW (SEE PARAGRAPH 6.2) 


Figure 6.29. RAVFAC Input With No Local Coordinate Systems Specified 


RAVFAC input also allows the user to specify an element relative to a local 
coordinate system. If one or more local coordinate systems are specified in 
the input, then one more card describing each one of these systems must be 
input. These cards relate each local coordinate system to the basic coordi- 
nate system. In this case the group of cards in RAVFAC input would be 
arranged as seen in Figure 6.30, as opposed to the manner shown in 
Figure 6.29 where no local coordinate systems were specified. 
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Figure 6. 30. RAVFAC Input With Local Coordinate Systems Specified 

After the last surface description card the user must insert a control card 
with -1 punched in columns 4 and 5 to identify the end of the surface descrip- 
tions. If the user describes one or more surfaces relative to a local coordinate 
system, then the program will expect a local coordinate system card for each 
different system referred to. After the last local coordinate system card, the 
user must insert another control card with -1 punched in columns 9 and 10. 
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Columns 


V ariable 
Name 



1 to 5 ISF Right- justified integer 5 0. First element number 

on the surface. Each surface may be divided into 
several elements. If there are two or more ele- 
ments on the surface, they will be numbered con- 
secutively from ISF. The numbering scheme can 
be seen in Figure 5. 8. Elements on different 
surfaces may be assigned the same number if 
desired. If ISF = 0 or is blank, the identification 
number will be one plus the identification number 
of the last element on the preceding surface. 


6 KSHD Integer 10. If = 0 the program is told that each 

element of the surface cannot cause shading. If 
> 0 each element of the surface is flagged as "can 
shade. " 

7 KBSHD Integer 10. If = 0 the program is told that each 

element of the surface cannot be shaded. If > 0 
each element of the surface is flagged as M can be 
shaded. " 

51 to 80 COMM Any comment can be made in these columns to identify 

the surface. 


6-31 



6,4.2 RAVFAC Card Number 2— Format 5X, 515, 5E10.5 
Variable 

Columns Name Description 

9, io ILK Right-justified integer indicating the surface shape. 


— ± 1 — 

Re c tangle 

= ±5 

— Cone 

= ±2 - 

Disk 

= ±6 

— Sphere 

= ±3 - 

Trapezoid 

= ±7 

— Circular 

-±4 — 

Cylinder 


paraboloid 


See section 5 for figures of each surface. For the 
flat surfaces (±1, ±2, ±3), the positive value indi- 
cates the +z side is active while the negative value 
flags the -z side as active. For the surfaces of 
revolution (±4, ±5, ±6, ±7), the positive value flags 
the outside of the surface as active while the nega- 
tive value flags the inside. If the user wants to 
specify both sides active he must input two separate 
surfaces, one looking each way. 

11 to 15 NVB Right- justified integer >0. Number of elements in 

the p direction. 

16 to 20 NVG Right- justified integer >0. Number of elements 

in the Y direction. 

21 to 25 NB Right- justified integer >0. Number of subelements 

per element in the p direction. 

26 to 30 NG Right-justified integer >0. Number of subelements 

per element in the Y direction. 

31 to 40 A Surface dimensioned (see figures in Section 5). 

41 to 50 BMIN Surface dimension Pmm ( see figures in Section 5). 

51 to 60 BMAX Surface dimension P max (see figures in Section 5). 

61 to 70 GMIN Surface dimension Y m i n (see figures in Section 5). 

71 to 80 GMAX Surface dimension Y max (see figures in Section 5). 
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6.4.3 RAVFAC Card Number 3 -Format 5X, 15. IPX, 6E10.5 


V ariable 

Columns Name Description 

6 to 10 NCS Right-justified integer. If NCS > 0 then the surface 

coordinate system is related to a local coordinate 
system numbered NCS. (Described on later cards. ) 
If NCS SO, then the surface coordinate system is 
related to the basic coordinate system (see Section 
3.2). 


21 to 30 RX 

31-40 RY 

41-50 RZ 

51 to GO PHI 


Distance of the origin of the surface coordinate sys- 
tem to the origin of the basic coordinate system (or 
local coordinate system if NCS > 0) measured along 
the x-axis of the basic coordinate system (or local 
coordinate system). 

Same as RX except that the distance is measured 
along the y-axis. 

Same as RX except that the distance is measured 
along the z-axis. 

Coordinate rotation angle, 4 1 » in degrees, for rotating 
the basic coordinate system (or the local coordinate 
system if NCS>0) into the surface coordinate system. 
Rotate Y toward X about Z. 


61 to 70 PSI Same as PHI but rotates X toward Z about Y (given 

in degrees). 

71 to 80 OMEGA Same as PHI but rotates Y toward Z about X (given 

in degrees). 


Note that the basic coordinate system (or local coordinate system) is rotated 
into the surface coordinate system. The rotations must be made in the order 
<}>then4' then u> (Figures 6.31 and 6.32). 
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Figure 6.31. Basic Coordinate System Rotation 



Figure 6.32. Rotations Order 
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6.4.4 Local Coordinate System Card— Format 5X, 15, 10X,6E10.5 


Variable 

Columns Name 


Description 


6 to 10 ICS 


Right- justified integer >0. Local coordinate sys- 
tem number. 


21 to 30 RXR 


31 to 40 RYR 


Distance of the origin of the local coordinate system 
to the origin of the basic coordinate system measured 
along the x-axis of the basic coordinate system. 

Same as RXR but measured along the y-axis. 


41 to 50 RZR 


Same as RXR but measured along the z-axis. 


51 to 60 PffiR 


61 to 70 PSIR 


Coordinate rotation angle, 4>, in degrees for rotating 
the basic coordinate system into the local coordinate 
system. Rotates Y toward X about Z. 

Same as PHIR but rotates X toward Z about Y. 


71 to 80 OMEGAR Same as PHIR but rotates Y toward Z about X. 


Rotation of the basic coordinate system into the local coordinate system must 
be carried out in the order 4>» w. 
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6. 5 360/95 JOB CONTROL CARDS 

The program delivery tape contains a copy of the source deck, a load module, 
and a copy of the job control language (JCL) necessary to run the program. 

At the Goddard Space Flight Center where the program was developed, the 
JCL has been placed in the 360/95 procedure library, so that only the JOB and 
EXEC control cards need appear before the data. 

For jobs where restart capability is not required, the EXEC card should be 
punched as 

//bEXECbVIEW, REGION=XXXK 

where b denotes a blank and XXX is a three-digit integer greater than 115, 
defining the requested core space. The amount of the request is discussed 
later in this section. 

For jobs where restart capability is desired, the EXEC card should be punched 
as 

//bEXECbVIEW, REGION-XXXK, TAP E S E R= Y Y Y Y YY , T A P E DSN=R EST A RT 

The XXX explained previously must be identical to the region request of the 
run which generated the tape. YYYYYY is the tape number of the restart tape 
that the user will create. (Usually a blank tape is submitted with the job. ) 

If punched output is desired, the following JCL card must be inserted after 
the EXEC card: 


//FT07F00 lbDDbDSN=& &C ARDS, SYSOUT=B 



6.5.1 Region Request 


The following formula may be used as a guide for estimating the region request 
made on the EXEC card for a given case: 

X = 218 X (Number of elements in case) + 115K 
Y = 48 X (Number of subelements in case) + 115 X (number of 
elements in case) + 115K 
Z = 55 X (Number of grid points in case) + 115K 
Request = max (X, Y, Z) 

If several cases are being run at the same time, the user must make the region 
request sufficiently large to handle only the worst case, that is, the one which 
requires the most core space. 

6.5.2 Time Estimation 

There are several factors which determine the time necessary to process a 
case. Among these are: 

• Method of integration used (see description of CASE control card 
columns 25 to 32). 

• Use of shading flags (see description of columns 41 to 48 of the CASE 
control card, and description of the NASTRAN-type input data card 
$VIEW, and RAVFAC-type input card 1). 

• Subelement mesh sizes specified for each element (the overall number 
of subelements determines the number of view-factor calculations). 

• Ordering of input data. The user can save time by ordering input 
so that the elements which cause the most shading appear first in 
the case data. Also, time will be saved if NASTRAN-type data con- 
tinuation cards appear directly after the cards which they continue. 

Since so many variables are involved in computing view factors, it is impossible 
to give the reader any accurate guidelines on estimating central processing 
unit (CPU) and input-output (IO) time requests. However, while the user is 
gaining experience, the use of the program's restart capability will ensure 
that runs which may time-out are not wasted effort. 
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A few examples of runs made at GSFC may serve as an indication of time 
requirements (refer to Table 6.1). 


Table 6. 1 

Sample Runs /Times 


Sub- 

Elements 

Shading in 
Problem ? 

CPU Seconds 
GSFC 360/95 
Slow Core 

Comments 

914 

YES 

994 


> ! 

Program asked to figure 
which elements shade 

969 

YES 

234 



1187 

YES 

781 


Program told which elements 

721 

YES 

38 


shade 

504 

YES 

44 J 



1448 

NO 

58 



1000 

NO 

62 


Program told that there is 

624 

NO 

12 


no shading 

480 

NO 

9 



1 150 

NO 

3 J 




6-38 




SECTION 7 


OUTPUT 



7. OUTPUT 


The first page of output contains the case number, date and initiation time 
of run, the problem title from the case title card, an echo of the case control 
card, and a description of its variables. If the output is produced from a 
restart case the first page will be as described, while the next page will echo 
the case control card submitted with the restart tape (see Section 9). 

A semisorted NASTRAN bulk data echo will begin on the next page if there is 
any NASTRAN-type formatted input. At its completion there will follow an 
input summary indicating the number of each of the card types (CORD1, CORD2, 
CHBDY, GRID and so on). Warning messages may also appear followed by 
a summary of any RAVFAC-type formatted data which was read. 

The computed view factors are printed out on an element-to- element basis, 
along with the areas of each element and the sum of all the view factors seeing 
each element. 

A summary of important program statistics is printed out. First, the user is 
told if the program was terminated normally or abnormally. In the latter case 
error messages appear. If the program terminated normally, the user is then 
informed of the CPU time used (in seconds) to process the input data, and 
the time used to compute view factors. Then the user is informed as to how 
many view' factor computations were made by the contour integral techniques 
and how many were made by the finite difference technique. The user is next 
informed as to the number of elements, subelements, and grid points that the 
data generated, along with the maximum number of elements, subelements, 
and grid points that the region request will allow. Finally, the user is in- 
formed as to the number of bytes of core not used because of an excessive 
region request. Hence, if the problem should have to be rerun, or if another 
problem of similar size must be run, the region request can be reduced to the 
minimal necessary size. 

The user can call for the program to produce punch output in the form of 
RADLST and RADMTX cards. The manner of producing this output is de- 
scribed in paragraph 6. 5. 

7.1 RADLST OUTPUT DATA CARD 

This card contains a list of element identification numbers given in the same 
order as the columns of the RADMTX matrix. Figure 7. 1 is the format and 
an example of a RADLST card. 
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5 6 7 8 


RADLST 

E1D1 

EID2 

EID3 

EID4 

E1D5 

EID6 

EID7 

EID8 

abc 

RADLST 

10 

20 

30 

50 

31 

41 

51 

61 

-hi oooooo 


+bc 

E1D9 

etc. 








+1000000 

71 










Figure 7. 1. Format and Example RADLST Card 
The contents of each RADLST output data card are as follows: 

Field Contents 

RADLST The word RADLST is punched beginning in column 1. 

EiDi The element identification numbers of the elements, 

given in the order that they appear in the RADMTX 
matrix (paragraph 7, 2) 

abc Sequential continuation numbers are punched out 

beginning with +1000000. 

7. 2 RADMTX OUTPUT DATA CARD 

This card contains numbers, described in Figure 7.2 and the following 
paragraph, so as to form one column of a matrix of radiation exchange co- 
efficients. Figure 7. 2 is the format and an example of a RADMTX card. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RADMTX 

INDEX 

R i j 

R i+1,i 

R i+2,i 

R i+3,i 

R i+4,i 

R i+5,i 

R i+6,i 

abc 

RADMTX 

3 

0 . 

9.3 

17.2 

16.1 

A 

0 . 

6.2 

+5000000 


+bc 

R i+7,i 

etc. 








+5000000 

6.2 










Figure 7. 2. Format and example RADMTX Card 


7-2 































The contents of each RADMTX output data card are as follows: 


Field 


Contents 


RADMTX The word RADMTX is punched beginning in column 1. 

INDEX The column number i of the matrix (Integer > 0) . 

R i+k t The matrix values (Real), starting on the diagonal, 

’ continuing down the column. K = 0, 1, 2, . . . , n - i. 

These values are the view factors from element i + k 
to element i multiplied by the area of element i + k. 

abc Sequential continuation numbers are punched out 

beginning with +5000000. 


The INDEX numbers go from 1 through NA, where NA is the number of radi- 
ating elements. Since the radiation exchange coefficient matrix is symmetric, 
only the lower triangle is output. Column 1 is associated with the element first 
listed on the RADLST card. Column 2 is for the next, and so on. 

NA 

P r = E R ij 0j 

i = i 

A = total irradiation at i 


qj - radiosity per unit area at j 

Rj j = view factor from element i to element j multiplied by the 
area of element i. 


The following three RADMTX and one RADLST card would form the lower 
triangular portion of a symmetric matrix as shown: 


RADMTX 

1 

0.0 

8.1 

RADMTX 

2 

3.2 

1.7 

RADMTX 

3 

12.2 


RADLST 

10 

50 

33 


"0.0 

0 

0 

R = 

8.1 

3.2 

0 


^ 2,2 

1.7 

12.2 


The first column describes the radiation exchange coefficients of element 10. 
The second column describes the radiation exchange coefficients of element 
50, and the third column, the radiation exchange coefficient of element 33. 
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8. ERROR AND WARNING MESSAGES 


Two types of messages may be printed out by the program at any point of the 
output. These are fatal errors and warning messages. A fatal error message 
terminates execution of the case being run. The program will usually try to 
run subsequent cases if they are present in the data deck. Warning messages 
apprise the user of a possible error in the data, but the program, being un- 
sure of the user's intent, will continue processing. 

In all cases, it has been attempted to make the messages sufficiently clear to 
explain the problem which caused them. However, since the explanations may 
be clear to the programmer, but may not be clear to the user, this section will 
expand on these explanations and/or refer the reader to the appropriate por- 
tions of this document. 

8. 1 EXPLANATIONS OF ERROR MESSAGES 

Fatal 

Error 

Number Explanations 

1 See Section 4 and paragraph 6. 5. 

2 The CHBDY card defines a line element. See Figure 6.8 and 
associated comments. 

3 Either the $VIEW card number in field 9 of the CHBDY card is mis- 
punched, or the $VIEW card is mispunched (field lor 2) or misplaced. 

5 Elements must be given positive areas. 

6 Coincident grid points will not properly define the element shape 
specified in field 4 of the CHBDY card. See Figures 6. 7 through 
6. 11 and associated comments. 

7 See Figure 6. 11 and associated comments. 

8 Fields 6, 7, and 8 of the continuation card of the CHBDY card were 
all found to be zero. Thus, a vector V = (0, 0, 0) was defined. 

The element type (either POINT or LINE) defined in field 4 of the 
CHBDY card requires V to be nonzero. See Figures 6.7 and 6.8 
and associated comments. 
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Fatal 

Error 

Number 


Explanations 


9 Either the GRID card number on the CHBDY card is mispunched 

or the GRID card is mispunched (field 1 or 2) or misplaced. 

10 All four grid points of AREA4 element must be coplanar. See 
Figure 6.11 and associated comments. 

11 The first three grid points, having been found collinear, do not 
properly define an AREA3 or AREA4 element. See Figures 6. 10 
and 6.11 and associated comments. 

12 The element type defined in the CHBDY card (either POINT or 
LINE) must reference a PHBDY card. The PHBDY card may have 
been mispunched or misplaced, or possibly the number in field 3 
of the CHBDY card was mispunched or missing. 

13 The CHBDY card does not define a REV element properly. See 
Figure 6.9 and associated comments. 

14 The CHBDY card defines either a POINT or LINE element. To 
either element, a vector V must be associated. The vector V 

is in the coordinate system defined in field 7 of the GRID card G1 
(specified in field 5 of the CHBDY card). That coordinate system 
was not defined by some type of CORD card. It was either missing 
or mispunched in fields 1 or 2. 

15 CHBDY cards defining POINT or LINE elements must have PHBDY 
cards associated to them. See description of CHBDY card, 

field 3, paragraph 6.3.1. 

21 Field 11 of a necessary continuation card was either misplaced or 
mispunched, or possibly field 10 of the card to be continued was 
mispunched. 

22 Field 3 of the CORD2 card must be blank! See description of 
CORD2 cards, paragraphs 6.3.5 through 6.3.7. 

23 User does not have the minimum amount of data necessary to de- 
fine one element by NASTRAN-type input data. Perhaps the 
ENDDATA card was misplaced. 
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Fatal 

Error 

Number Explanations 

24 There must be an ENDDATA card if the case has NASTRAN-type 
input data. See paragraph 6. 1. No subsequent cases are processed. 

25 Too many elements have been defined and core space has run out. 
Increase the region request or lower the number of elements in 
the model. See paragraph 6. 5. 

26 More than 40 CORD cards of either type are not permitted. 

27 Field 4 was probably mispunched. 

28 In field 1 of the CORD card, CORDXY must be punched where 

X = 1 or 2 and Y = R, S, or C. The card had X or Y mispunched. 

29 Field 11 of each continuation card must be unique. The deck has 
at least two with the same number. 

40 Either field 3 of the GRID card has been mispunched or the CORD 

card which was intended to be referenced has been mispunched 
(field 1 or 2) or misplaced. 

50 

51 Probable mispunch. 

52 

53 

54 No embedded blanks are allowed in data. 


> Should not occur. Seek a programmer’s assistance. 

56 ) 

60 Check the region request on the generating run. Use the same 

request on the restart deck JCL. Resubmit. 

611 l May have a bad tape. Try once more. 

621 l If not successful, start over again. 

70 The grid points used to define coordinate systems by means of 

CORD1 cards must be in the basic coordinate system. 
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Fatal 

Error 

Number Explanations 

80 Enough space was requested to load the program, but not enough 
to store any data. Increase the request. See paragraph 6. 5. 

81 Error should not occur. If it does, there is an error in subroutine 
GETMA. You will need a programmer's assistance. Save the 
run to show the programmer. 

90 The RAVFAC surface shape that was specified in columns 9 and 10 

of a RAVFAC number 2 card is undefined. It must be either ±1, 
±2, ±3, ±4, ±5, ±6, ±7. 

93 The region request was not large enough for all the elements re- 

quested, Increase the region request or decrease element usage. 
See paragraph 6.5. 

100 See paragraph 6. 1. 

101 The region request was not large enough for all the elements or 
subelements requested. See paragraph 6.5. 

102 \ 

110 > See paragraph 6. 1. 

111 ) 

112 See description of TITLE card, paragraph 6.2. 

113 See description of CASE card, paragraph 6. 2. 

114 The CPU time request must be increased. Consider using restart 
capability next time. See paragraph 6.5. 

8. 2 EXPLANATIONS OF WARNING MESSAGES 


Warning 

Number Explanations 

4 The VIEW card has a blank or zero in fields 5 and/or 6. This does 

not properly define a mesh of subelements so the program assumes 
a one-by-one mesh (NB = 1 and NG = 1) and continues processing. 
Check that this did not affect the answers. 
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Warning 

Number 


Explanations 


16 See paragraph 2.5. 

30 The program gets tired after checking 500 out-of-order continua- 
tion cards for duplicate ID fields. It refuses to check any remain- 
ing continuation cards. The user should check. 

31 See description of GRDSET card in paragraph 6.3. 9. 

32 If continuation cards do not appear directly behind the card they 
continue, then processing is slowed. 

33 Some card types were mixed in the NASTRAN-type input data 
which were not used by VIEW. It is wise to check that these are 
not mispunched cards which were really intended for VIEW 
processing. 

34 The NASTRAN-type input data element ELCYL is not recognized, 
by VIEW. Other element types will have to be used to approximate 
the ELCYL element. 

35 In the CHBDY card, which specified a POINT or LINE element, 
the continuation field (10) was left blank. The program interprets 
this as implying that the user does not want this element to appear 
in his model for calculation of view factors. 

91 The RAVFAC-type input card number 2 has a blank or zero in 
columns 11 to 15 and/or 16 to 20. This does not properly define 
a mesh of elements so the program assumes a one-by-one mesh 
(NVB and NVG = 1). Check that this did not affect the answers. 

92 The RAVFAC-type input data card number 2 has a blank or zero 
in columns 21 to 25 and/or columns 26 to 30. This does not 
properly define a mesh of subelements so the program assumes a 
one-by-one mesh (NB and NG = 1) . Check that this did not affect 
the answers. 

103 No further cases are run. 
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9. RESTART 


Restart is a program feature which saves completed work on a tape; should 
disaster occur, the user can resubmit his restart JCL with the restart tape 
and pick up where he left off. 

For example, suppose a user has just put together a large model which he ex- 
pects will take no longer than 30 CPU minutes and 10 I/O minutes to run. As- 
sume that one of the following happens: 

• The job really takes 31 CPU minutes to run. 

• The job really takes 11 I/O minutes to run. 

• The computer operator accidentally hits the emergency OFF switch. 

Normally such events mean that whatever calculations have already been done 
are wasted because the whole job will have to be run over again. Restart 
solves this problem. 

To create a restart tape during the original run the user should: 

• Punch any positive integer anywhere in columms 17 to 24 of the case 
control card (see paragraph 6.2). 

• Have the JCL EXEC card read 

//bEXECbVIEW, REGION=XXXK, TAPESER =YYYYYY, TAPEDSN=RESTART 

where b denotes a blank and XXX should be the region request and 
YYYYYY is the tape number (see paragraph 6.5). 

To make a run from a restart tape the user should: 

• Have the same EXEC card just described. 

• Have as input data a single card with RESTART punched in columns 

1 to 7, and anywhere in columns 9 to 16 an integer which is the same 
number as the CPU time estimate on the JOB card. If the time 
estimate on the job card is one-half minute, punch an H anywhere in 
columns 9 to 16. 
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Submit the JCL, single data card, and the restart tape. 


The data for more cases may be placed behind the previously described single 
data card. 
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APPENDIX A 


SAMPLE PROBLEM USING ONLY NASTRAN-TYPE DATA 


The physical structure to be modeled in this problem is shown in Figure A-l. 



Figure A-l, Structure to Be Modeled 

As shown in Figure A-l, nonzero view factors exist only between the two 
sides of plate A and plate B which face each other (call them face A and face B), 
and between the outside of the cylinder (face C) and faces A and B, Hence, the 
model, if desired, could consist of only five elements: two rectangular, one 
cylindrical, and two disks. The disks would cap the ends of the cylinder so 
that the two rectangles would not see each other through the cylinder. 

Choosing the axis of the cylinder as the z-axis and the larger rectangle to be 
in the x-y plane, the element model in the basic coordinate system can be 
drawn as seen in Figure A-2. 


A-l 



ELEMENT 30 



Figure A-2, Element Model in the Basic Coordinate System 

The grid points are arbitrarily numbered as shown for the two rectangular 
plates. Figures A-3 and A~4 show more clearly the grid points used to define 
the cylinder and two disks. 



Figure A-3. Grid Points Used to Figure A-4. Grid Points Used to 

Define Cylinder Define Two Disks 


A-2 


Since shading is involved, a fine subelement mesh is needed on elements 10 
and 30, and a less fine subelement mesh is needed on element 20. Since ele- 
ments 40 and 50 are present for shading purposes only, we can specify a one- 
by-one mesh on each of these elements. Furthermore, we can specify that 
element 10 can be shaded but cannot cause shading. The same is true for 
element 30. However, element 20 can shade but cannot be shaded while 40 
and 50 can shade and be shaded. 

Figure A-5 is a listing of the input used to run this problem. It is followed by 
a copy of the program output. 


//B9EFPT01 JOB (W1831C* T »B0008I » 00 1H00 > * 1 YY »MSGLEVtL= <2*0 > 

// EXEC VIEW,REGI0N=175K 

TITLE SAMPLE PROBLEM NO. I 

CASE= 1 

BEGIN BULK 


CHBOY 

10 


AREA4 

10 

11 

12 

13 

50 

CHBDY 

20 


REV 

20 

30 



55 

CHBDY 

36 


AREA4 

100 

103 

102 

101 

60 

CHBOY 

40 

10 

POINT 

50 



65 

♦ Cl 

♦ Cl 





0.0 

0.0 

1 . 


CHBOY 

50 

10 

POINT 

40 



65 

♦C2 

♦C2 





0.0 

o 

. 

o 

-1.0 


SVIEW 

50 


0 1 

25 

15 




SVIEW 

55 


1 0 

10 

20 




SVIEW 

60 


0 1 

15 

15 




SVIEW 

65 


1 1 

5 

5 




PHBOY 

10 


.786 






GRID 

10 


2.5 

-5.0 

0.0 




GRID 

11 


2.5 

5.0 

0.0 




GRID 

12 


-2.5 

5.0 

0.0 




GRID 

13 


-2.5 

-5.0 

0.0 




GRIO 

100 


2.5 

-2.5 

10.0 




GRID 

101 


2.5 

2.5 

10.0 




GRID 

102 


-2.5 

2.5 

10.0 




GRIO 

103 


-2.5 

-2.5 

10.0 




GRID 

20 


.5 






GRID 

30 


.5 


10.0 




GRID 

40 








GRID 

50 




10.0 





ENOOATA 

ENDCASE 


Figure A-5. Program Input for Sample Problem 1 
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A-4 


G 0 0 0 A R 


THE 



ED PUCCINELLI 
REG MITCHELL 
CLIFF JACKSON 


VERSION ONE OCTOBER 31* 1972 


Q * 


FLIGHT CENTER NASA 

»»»§»« 4 « 4 « 


EEEEEEE 

E 

£ 

E 

EEEE 

EEEE 

E 

E 

E 

EEEEEEE 



W Ww 


PROGRAM 
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DATE 

TIME 




* CASE l * 


PROBLEM TITLE *** SAMPLE PROBLEM NO. 1 *** 


CASE CONTROL CARD 


COLUMNS 

VARIABLE 

VALUE 

DESCRIPTION 

9-16 

IENDTM 


1 

ENTER CPU TIME USED ON J06 CARD <IN MINUTES— USE H FOR ONE-HALF MINUTE) 
DEFAULT IENUTM=iOOOOO (SEE DOCUMENTATION). 

17-24 

NT 


0 

LE 0 
GT 0 

-DO NOT CORY INPUT DATA ONTO TAPE (UNIT 2) FOR RESTART USE 
-COPY INPUT DATA ONTO TAPE (UNIT 2) FOR RESTART USE 

25-32 

NVFCAL 


0 

GT 0 
LT 0 
EQ 0 

-VF CALCULATED USING FINITE DIFFERENCE METHOO 
-VF CALCULATED USING CONTOUR INTEGRATION METHOD 
-PROGRAM SELECTS METHOD TO BE USED (SEE DOCUMENTATION) 

33-40 

NFE 


0 

LE 0 
GT 0 

-MESH SIZE FOR EACH ELEMENT SPECIFIED ON INPUT DATA 

-MESH SIZE SET TO 1 BY i FOR EACH ELEMENT. (OVERRIDES INPUT DATA) 

41-48 

NFS 


0 

GE 0 
LT 0 

-SHADING CONSIDERED USING DATA FROM INPUT 
-NEGLECT ALL SHADING 

49-56 

RMAX 

0.0 


MAXIMUM AREA/OISTANCE RATIO. OEFAULT RMAX * 0.1 (SEE DOCUMENTATION). 


7 / 26/73 

14 . 44.40 
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nastran formatted input 


SEMl-SORTEO BULK DATA ECHO 


1 

2 

3 


<* 

5 



«««««««» 



CHBDY 

10 



AREA4 

10 

CHBDY 

20 



REV 

20 

CHBDY 

30 



AREA4 

100 

CHBOY 

40 

10 


POINT 

50 

♦ Cl 






CHBDY 

50 

10 


POINT 

40 

♦ C2 






GRID 

10 



2.5 

-5.0 

GRID 

11 



2.5 

5.0 

GRID 

12 



-2.5 

5.0 

GRID 

13 



-2.5 

-5.0 

GRID 

100 



2.b 

-2.5 

GRID 

101 



2.3 

2.5 

GRID 

102 



-2.5 

2.5 

GRID 

103 



-2.5 

-2.5 

GRID 

20 



• 5 


GRID 

30 



.5 


GRID 

40 





GRID 

50 





PHBOY 

10 



.786 


SVIEW 

50 


0 

1 

25 

$V 1 EW 

55 


1 

0 

10 

SV1E# 

60 


0 

1 

15 

SVIEW 

65 


1 

1 

5 


6 

7 

8 

9 

10 




««»«&««« 

««««««»« 

11 

12 

13 

50 


30 



55 


103 

102 

101 

60 





65 

♦Cl 

0.0 

0.0 

1. 

65 

♦C2 

0.0 

0.0 

-i.o 




0*0 

0*0 

0*0 

0*0 

10*0 

10.0 

10.0 

10.0 

10.0 

10.0 


15 

20 

15 

5 


CONTINUED 
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NASTRAN FORMATTED INPUT 
SEMI-SORTED BULK DATA ECHO 


INPUT SUMMART 

NUMBERS OF COMPLETE LOGICAL CARD TYPES DETECTED 

CH8DY = 5 CORD1 = 0 C0RD2 - 0 GRDSET = 0 GRID - 12 PHBDY * 1 VIEtf * A 



A-8 


ELEMENT 

10 

20 

30 

40 

SO 


VIEW FACTORS 


AREA 


ELEM TO ELEM = VIEW FACTOR 


ELEM TO ELEM = VIEW FACTOR 


VF SUM FROM 
ELEMENT 


0 « 50000E Q2*******#*************************************************************************** 1 **** 
10-20 =0 .673263E-0 1 20-10 -6.107153E 00 

10 - 30 *0 *5Q3031E-0 1 30 - 10 =Q.10G606E 00 

0 » 1 1763E 00 

20 - 30 =0.728946£-Ql 30 - 20 =6 .916021E-0 1 

0,ltt005E 00 

■ “ ■ ... 0 * 19221E 00 

0* 7860 OE Q0***** tt *#*** # #****************’ t ***** <> **** fl ' < ******^*' > * ,,>o ' &<> ******* tt ******* ******************** 

o.o 


0*7860^E 00 * **«~***f>* 


■i.'* 
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PROGRAM SUMMARY FOR CASE l 


PROGRAM TERMINATED NORMALLY 


CPU TIME (SECONDS! 


PROCESSING OF INPUT DATA- 0 

VIEW FACTOR COMPUTATIONS* 19 


NUMBER OF VIEW FACTOR COMPUTATIONS 


FINITE DIFFERENCE* 115221 

CONTOUR INTEGRAL* 590 


REQUESTED CORE SPACE ALLOWS THE INPUT DATA GENERATED 


300 

ELEMENTS 

5 

ELEMENTS 

660 

SUB-ELEMENTS 

850 

SUB-ELEMENTS 

1360 

SUB-ELEMENTS 




AFTER COMPRESSION 



900 

GRIO POINTS 

12 

GRID PO IN ) S 


21K BYTES OF CORE WERE NOT USED IN THIS CASE. 


( 62K ON ELEMENTS* 21K ON SUB-ELEMENTS AND 62* ON GRIO POINTS.) 


««««««* ft 

* END CASE 1 * 



APPENDIX B 


SAMPLE PROBLEM USING NASTRAN AND RAV FAC-TYPE DATA 
The physical structure to be modeled in this problem is shown in Figure B-l. 



Figure B-l might be an approximation of a spacecraft in a thermal vacuum 
chamber. Suppose it is desired to determine the view factors between the 
exterior of the spacecraft and the interior of the chamber. 

The model should be made using NASTRAN- type data for the spacecraft and 
RAVFAC-type data for the chamber. Locating the origin of the basic coordi- 
nate system at the center of the spacecraft, the element model in the basic 
coordinate system can be drawn as seen in Figure B-2. 


B-l 



z 



•5 


Figure B-2, E lenient Model in the Basic Coordinate System 
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Table B-l lists the optional shading flags. 

Table B-l 
Shading Flags 


Element 

Can Shade 

Can Be Shaded 

10 

YES 

NO 

11 

YES 

NO 

12 

YES 

NO 

13 

YES 

NO 

14 

YES 

NO 

15 

YES 

NO 

100 

NO 

YES 

101 

NO 

YES 

102 

NO 

YES 


Figure B-3 is a listing of the input used to run this problem. It is followed by 
a copy of the program output. 


//B9EFPT02 JOB <®B«lB31C*TtBO00ai*O09OO2> ♦YYY t MSGLEVtL= <2i 0) 

// EXEC VI£W*R£GI0N"200K 

TITLE /// SAMPLE PROBLEM NO. 2 /// 

CASE" 9 

BEGIN BULK 


GRID 


1 


1.0 

-1 

.0 

-1.0 





GRID 


2 


1.0 


1*0 

-1.0 





GRID 


3 


■*1.0 


1.0 

-1.0 





GRID 


4 


"1.0 

-1 

.0 

-1.0 





GRID 


5 


1.0 


1.0 

1.0 





GRIO 


6 


-1.0 


1.0 

1.0 





GRID 


7 


-1.0 

-i 

.0 

1.0 





GRID 


8 


1.0 

-1 

.0 

1.0 





CH0OY 


10 


AREA4 


1 

2 

5 

8 

25 


CHBOY 


11 


AREA4 


2 

3 

6 

5 

25 


CHBOY 


12 


AREA4 


3 

4 

7 

6 

25 


CHBOY 


13 


AREA4 


8 

7 

4 

1 

25 


CHBDY 


14 


AREA4 


4 

3 

2 

1 

25 


CHBOY 


15 


AREA4 


8 

5 

6 

7 

25 


SVIEW 
endoata 
RAVFAC 1 

DATA 

25 


1 0 


5 

5 





100 l 








END OF 

CHAMBER 



-2 

1 

1 

5 5 

7.5 


0.0 

b.O 

0.0 


360. 











-90. 

102 1 








OTHER 

END OF 

CHAMBER 


-2. 

1 

1 

5 5 

7.5 


0.0 

b.O 

0.0 


360. 












90. 

101 1 








chamber sides 




-4 

1 

1 

40 20 

5.0 


-7.5 

'.5 

0.0 


360. 


90. 

-I 

ENDCASE 


Figure B-3. Program Input for Sample Problem 2 
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B-4 


GOODARO SPACE FLIGHT CENTER NASA 
#******_*##•# *****♦.*_ * *__ . * *_* * 


THE 


7 

V 

V 

V 


V 

V 

V 

V 


V V 

V V 

V V 

V V 


vv 

vv 


EEEEEEE 

E 

£ 

E 

EEEE 

EEEE 

E 

E 

t 

EEEEEEE 


W W 

w W 

w. w 

W W 

M W M 
WWW 
W WWW w 
W WWW w 
WW WW 
WW WW 


PROGRAM 


EO PUCCINEllI 
REG MITCHELL 
CLIFF JACKSON 


VERSION ONE OCTOBER 31* 1972 



B-5 


QAIE 

TIME 


* CASE 1 * 


PROBLEM TITLE 

/// 

SAMPLE PROBLEM NO. 

2 /// 

CASE CONTROL 

CARO 



• — • - 

COLUMNS 

VARIABLE 

VALUE 

DESCRIPTION 

9-16 

IENDTM 

6 

ENTER CPU TIME USED ON JOB CARD (IN MINUTES— USE H FOR ONE-HALF MINUTE) 
DEFAULT IENDTM=iOOOOO (SEE. DOCUMENT AT ION).. ... _ ... -• 

17-24 

NT 

0 

LE 0 
GT 0 

-DO NOT COPY INPUT DATA ONTO TAPE (UNIT 2). FOR RESTART USE 
-COPY INPUT DATA ONTO TAPE (UNIT 2) fOR RESTART USE 

25-32 

NVFCAL 

0 

GT 0 
..LI 0 
EQ 0 

-VF CALCULATED USING FINITE DIFFERENCE METHOD 
. -VF calculated using contour ...integration method 

-PROGRAM SELECTS METHOD TO BE USED <SEE DOCUMENTATION) 

33-40 

NFE 

0 

LE 0 
GT 0 

-MESH SIZE FOR EACH ELEMENT SPECIFIED ON INPUT DATA 

-MESH SIZfc SET TO 1. BY 1 FOR EACH. ELEMENT . (OVERRIDES. INPUT. DATA). 

41-46 

NFS._ 

Q._._ 

. G£ 0 
LT 0 

-SHADING CONSIDERED. .USING DATA FROM INPUT 
-NEGLECT ALL SHADING 

49-56 

RMAX 

0,0 

MAXIMUM AREA/OISTANCE RATIO. DEFAULT RMAX = 0.1 (SEE OOCUMENT AT ION) . 


7/2T/J3. 

9.27.29 
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NAS I RAN PORKATTEO INPUT 
SENS -SORTED OULK DATA ECHO 


1 

2 

3 . 4 . - 


5- 

- & 

-7 

a - 

9 

10 - - - 







' . 




— 

CHBDY 

10 

AREA4 


1 

2 

s 

6 

25 


CHBDY 

1 1 

. .. AREA4 

. _ . 

2. 

3 

.. . .6. 

5_. . - 

25 

— 

CHBDY 

12 

AREA4 


3 

4 

7 

6 

25 


CHBDY 

13 

. . AREA4 „ „ 


8 

.7 

4 

1 

25 

- • 

CHBDY 

14 

AREA4 


4 

3 

2 

1 

25 


tHBDY . . 

15 

AREA 4 * ... 


8 

5 

_6 _ ... 

7 

25 

*■• • 

GRID . 

1 

l.Q 


-1*0 

-1.0 ... 

— 





GRID 

2 

U 

0 

1.0 

-UO 





GRID 

3 

-1. 

0. 

1.0 

^UQ 

- • 

— 



GRID 

4 

-uo 


-1.0 

-1*0 





GRID 

5. 

. ... .1a. 

0 

. 1*0_. 

. _ 1*0 

— - — - 


- . 


GRID 

6 

-1. 

0 

uo 

uo 





GRID 

7 . 

-1*0 


-1.0 

. uo ... 

„ .... — 




GRID 

8 

1.0 


-1*0 

uo 





$view 

_ ' 25 

1 0 


5 

5 








_ 


.. .. ... . 


- - - • • - 

. ... .. 

— - 





... INPUT _£UMKARV 


CH60Y = 


Ji UMBERS OF COMPLETE LOGICAL. CARCL.TYPES J1EIGPTE0-. - - 
CORD.i - 0„_. C0RD2. = ... 0. 3RD SET .= . ,0 OR I D — 8. 


PHBDY 


__VI£¥.=. 





RAVFAC DATA 


SURFACE INPUT DATA 


SURFACE NUMBER 100 CAN ShADE - NO END OF CHAMBER 

... - CAN BE- -SHADED- YES . _ . 

SURFACE TYPE * -2_ .ALPHA .? 0.7504 OE 01 

NODES* B-DIR.= 1 ELEM/NODE. B-DIR.= 5 BETA(MIN) = 0.0 BETA(MAX) = 0.50000E 01 

NODES. G=OIR*=« 1 ELEM/NODE* G-01R,=. S GAMMA! MINI .=... 0 .Q . GAMMA (MAX) =_ _fl_, 36QO0E 03 - 

RXXL5_(L.iL R(Y1=. 0^0 iliZl = 0.0 ... 

PHI * 0.0 PSI = 0.0 OMEGA=-0 . 90000E 02 REF. TO INTERMEDIATE COORDINATE SYSTEM 0 


w 

I 

-O 


SURFACE NUMBER 200 CAN SHADE - NO OTHER END OF CHAMBER 

CANBESMADED- YES ' _ ... 

SURFACE TYPE = . -.2 ALPHA = D.2S00QE.D.1 

NODES* B-OIR.= 1 ELEM/NODE. B-DIR.= 5 BETA (MIN) = 0.0 BET A (MAX ) = 0.50000E 01 

NODES G-OIR* = ». 1 ELEM/N0QF*_G=D1R.= S_ _ ... £AMMA.(M1NJ = 0.0 . GAMMA ( MAX 1= 0.3A00QF 01 ..... 

RXX)=0.0 — iUYI= 0.0 J?.(Z1. =.0.0 ... 

PHI = 0.0 PSI = 0.0 0MEGA= 0.90000E 02 REF. TO INTERMEDIATE COORDINATE SYSTEM 0 


SURFACE NUMBER 300 CAN SHADE - NO CHAMBER SIDES 

... . CAN. BE .SHADED- _ YES 


SURFACE TYPE = -A . ALPHA. = 0.S0000E 01. . . 

NOOES. B-DIR.= l ELEM/NODE* B-DIR.= AO BET A (MIN) =-0.7500QE 01 BETA(MAX) = 0.75000E 01 

NODES. G-DIR.S* 1 ELEM/NOQ£*_G-D I R. _=_... 20. GAMMA (MIN I » 0.0 GAMMA ( MAXI* 0. 36000g 03 

JUX) «. 0.0 ... ... R ( Y) - 0 .0 fiXZJ = .0.0 _ 

PHI = 0.0 PSI = 0.0 OMEGA* 0.90000E 02 REF. TO INTERMEDIATE COORDINATE SYSTEM 


0 
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VIEW FACTORS 


ELEMENT 


AREA 


ELEM TO ELEM = VIEW FACTOR 


ELEM TO ELEM 

_«.« *. • • • ■ ^ m 


VIEW FACTOR 


VF SUR E-SQfrL 

ELEMENT 


LM-iJJ *.*-■•»* - 


i n ft.anQQDF ft l »»«»»»» tt»»tt »tt»»tt»»»»»»«» »»*»»*»***** ***»********** ***** *********** ******^^P*-*~-^^— - 

10 - 100 *0.2?245uE-01 100 - 10 =0* 138757E-02 

14. • 200.. *G^72*4SEr_0 1 .20<L-.__10 =0.13fi75IE-02 .... 

10 - 300 =0 « 944122E 00 300 - 10 *6 .80 1395E-02 

0*9986 IE 00 


11 Q.40QQ0E 01*******^*********** «*#*»***»** ********************************************************* 

" 11 - 100 =0.368j80E 00 100 - 11 =6 . 187614E-0 1 

11 - 300 =0.fi326(lM: .00 300...... .11. =fi.»53t971E-02 --- 

■ " ■' 0*100 10E 01 

12 0* 4Q000E 0 1 ************************** ************************************************************* 

# _ „ _ iz - UOQ =Q.^72449ErflL .._. 1Q0 - _! _ 12. =0^1 38157 Er 02 

12 - 200 =0.272449E-01 200 - 12 =0 * 1 38757E-02 

_ 12 - _3Q.O =0.944l22E .00 ... 300 - 12„ ?0 -3Q 1395E-02 


0.99661E 00 


13 0.4QOOOE 01**********************^********* ^t***-*********** »******>****** »*&***« 

13- 200 =O.J68iOOL QO ... . . 20 j 5 -' _.l3 =0.187.61' ^-01. . 

300 - 17 *0.536971E-02 

_ _ _ . j).iooij)e_oj_. 


13 


300 =0.fo32b04E GO 


Q ■ 4.0 D Q 0£ 0.1 *. ****** ****^**^* ********!.******** * * * * ** * *** * ** »»*»» » a » »#» ft ® < *».** ******* *******L 

“ " 14 - 100 =0.29i990E-01 100 - 1^ =S.16870‘3E-02 

14 - 200 aQ.233474E-01 200 _r... . 14. _=Q .U.0.90.7 E“0 2 

14 - 300 =0.944i22E 00 300 - 14 =0.801395E“02 

.0.9.96 ME JWL- 

^5 q ^qoOOlE 01 ^-fr****** *************** ********** ********** *********_' **J****** # ** *********** ****' *********_ 
*---■— — fs^" ioo =0^2334 ?5E-01 ' 100 - IS =0. 118908E-02 

15_r 200 =Q.,291?$L1E-01 20 0_- . . 15 . «0-*1*&7.10£.-02. 

15 - 300 =0 .9441 14E 00 300 - 15 =0.301388E-02 

__ . 0.9 9666 E M.. 

1 Art A 7QCAAC fl?***************************************************************************^******^-***^ 

’ Too - — 2‘00 : =0.“766il7E-0r 200‘- 100 =0.7$oVl6E-of 

100^ 300_=fi,892820£_00 J3 M_t .. .10.0. ”0^.1468.0 3.E. .00. 

- -• - - 0.99204E 00 


200 0.78540E 


_2DO_.T- 31L0_^ft» 892821£_0O — 1 


366 - 200 . = (L.1 GBB04E 66 ' 


0.99205E 00 


0.99579E 00 
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PROGRAM SUMMARY^ FOR- CASEU 

«»»***««*»»««»««»•«•*«*•»« 


PROGRAM TERMINATED NORMALLY 


CPU TIME < SECONDS) 


PROCESSING OF. INPUT DATA*. 0 .. 

VIE* FACTOR COMPUTATION 166 


NUMBER OF VIEW FACTOR COMPUTATIONS 


FINITE DIFFERENCE 3 651256 
CONTOUR INTEGRAL 3 1390 Z.„ 


REQUESTED CORE SPACE ALLOWS 


THE INPUT DATA GENERATED 


416 ELEMENTS 9 ELEMENTS 

094 SUB-ELEMENTS 1000 SUB-ELEMENTS 

1862 SUB-ELEMENTS 

AFTER COMPRESSION 

1240 GRID POINTS 8 GRJO POINJS 


39K BYTES QF CORE WERE NOT USED IN THIS CASE, 


< „ 00K.ON ELEMENTS* 3VK ON SUB- ELEMENTS AND 8/K ON_.GRJD, POINTS , ) 


*««»««»«» «•««« 

* _ END ' CASE. \_* 



APPENDIX C 


SAMPLE PROBLEM DEMONSTRATING THE USE OF LOCAL 
COORDINATE SYSTEMS 


The physical structure to be modeled in this problem is shown in Figure C-l. 



2X2 SQUARE 
PLATE B 
PARALLEL TO 
SQUARE PLATE A 


Figure C-l. Structure to Be Modeled 

The structure to be modeled consists of the two opposite sides of a cube 
"looking" at each other. Plate A will be modeled using NASTRAN-type data 
and using a local coordinate system while plate B will be modeled using 
RAVFAC-type data and another local coordinate system. 

The basic coordinate system will be chosen to lie halfway between the two 
plates with the X-Z plane parallel to the plates. The local coordinate systems 
will be simple translations and rotations so that the X-Y planes contain their 
respective plates. This is shown in Figure C-2. 


C-l 



z 



Figure C-2. Element Model and Its Coordinate Systems 

The coordinates of grid points 10 through 13 with respect to local coordinate 
system A are respectively (1, 1, 0), (-1, 1, 0), (-1, -1, 0), (1, -1, 0). 

Figure C-3 is a listing of the input used to run this problem. It is followed by 
a copy of the program output. 


C-2 



//fl9£F PTOJ JOB (Ml*lb31CtT»B0008 WHOOHOOWYyWMSGL£Vt.L=<2,C> 
// EXEC VIEW »REG ION= 1 50K 

TITLE *** SAMPLE problem NO* 3 *** . . . .. . . . 

C.ASE= H 

BEGIN BULK 


CHBOY 

100 


AREA4 

10 

11 

12 13 SO 

_.SVIEW 

50 

0 

0 

2 

2 

.... 

GRID 

10 

s 

1.0 

1.0 

0.0 


GRID 

11 

s 

-1.0 

1.0 

0.0 


GRID 

12 

s 

-1.0 

-1.0 

0.0 


GRID.. 

13 

s 

1.0 

-1.0 

0.0 

0.0 -.5 WO *C1 

C0RD2R 

S 


0.0 

-1.0 

1.0 

_* CJL 

*5 

-.5 

1.0 




ENODATA 







BAVJEACL. 

data 





..... . . . . 

20 






PLATE B 


♦1 1 

1 2 

2 

0.0 

-wo 

1.0 -1.0 wo 


15 






=1._ 






90.0 


"is 

0.0 


1.0 1 

.0 


_=1 




_ 





ENDCASE 


Figure C-3. Program Input for Sample Problem 3 


C-3 



C-4 


GOODARU SPACE FLIGHT CENTER NASA 

ff » * W « « « **■«■*«■ 4 » tt * ft * 4 4 4 4#* **4* 


t nt 


v 

V 

V 

V 

V 

V 


1 / 

V 

V/ 

V 

V 

V 


V v 

V V 

vv 

vv 


EttEEEE 

E 

t 

E 

EEEE 

ElEE 

E 

t 

t 

EEEEEEE 


W v/ 

rf R 

M W 

W M 

w w w 

W rf tf 


* W*W W 

w ww# w 


W«J ww 
WW tow 


PROGRAM 


ED PUCCINELLI 
REG MITCHELL 
CLIFF JACKSON 


VERSION ONE UCTOdER 31* i'ild 
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TIME 1 1 *23*23 


««««»««»«« 
* CASE 1 * 


PROBLEM TITLE 

«<** SAMPLE 

PROBLEM MO. J »** 


CASE control 

CARD 






COLUMNS 

VARIABLE VALUE 

description 



9-16 

ienotm 

H 

ENTlR CPU time uSt D ON JOB CARD (IN MINUTES--USE H FOk ONE-HALF 
□EFAUlT ILNUTM=iUOOOO (StE DOCUMENTATION). 

' MINUTE f' 

17-24 

NT 

0 

LE 0 -DO NOT COPY INPUT DATA ONTO TAPE (UNIT 2) FOP WESTAHT USE 
GT 0 -COPY INPUT DATA ONTO TAPE (UNIT 2) FOR RESTART USE 


2S-32 

nvfcal 

0 

GT D -VF CALCULATED USING FINITE DIFFERENCE METHOD 

LT U - VF CALCULATED USING CONTOUR INTEGRATION METHOD 

EQ 0 -PROGRAM SELECTS METHOD TO BE USED (SEE DOCUMENT AT ION) 



33-40 

NFt 

0 

LE 0 -MESH SlZt FOR EACH ELEMENT SPECIFIED ON INPUT DATA 

GT 0 -MESH S IZl SET TO I BY 1 FOR EACH ELEMENT. (OVERRUNS INPUT 

DATAL 

41-4G 

NFS 

0 

GE U -SHAOlNG LUNSIDEREO USING DATA FROM INPUT 
LT 0 -NEGLECT ALL SHADING 



49-bb 

HMAa 0*0 


MAXIMUM AREA/DISTANCE RATIO. DEFAULT RMAX = 0.1 (SEE DOCUMENTATION). 
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N ASIRAN FORMATTED INPUT 
SEMI-SORTED BULK DATA ECHO 


1 

«««#«««« 

CMBDY 

2 

««««»«»« 

100 

3 

A 

ARLAA 

5 

«««««««« 

10 

6 .7 „ 9 ... 1Q_ 

».«.«««« «««•»«««« ««««»««» «»«««*»« 

11 12 13 50 

C0RD2R 
♦ Cl 

s 

.5 

-.5 

0.0 

1.0 

-1.0 

1.0 0.0 -.5 1.0 ♦Cl 

GRID 

10 

5 

1.0 

1.0 

0.0 

GRID 

11 

S 

-1.0 

1.0 

0.0 

GRID 

13 

s 

-1.0 

-1.0 

0.0 

GRID . 

13 

s 

1.0 

-1.0 

0.0 

SVIEW 

50 

0 

0 

2 

2 


INPUT SUMMARY _ __^ 

NUMBERS OF COMPLETE LOGICAL CARO TYPES_J)ETECTEO 

CHBDY - 1 CORD 1 - . 0 C0RD2 = 1 GROSET = 0 . G R ID = A P HBDY = 0 VIEW = 1 
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RAVFAC DATA 


SURFACE INPUT DATA 


SURFACE NUMBER 

20 

CAN 

SHADE 

NO 

PLATE B 



CAN 

BE SHADED- 

NO 




SURF ACE TYPE = I 

•NODES* b-OIR.= 1 ELEM/NODE* B-DIK.= 

NODES G“D1R*=* 1 ELEM/NODE* G-DI*.= 


ALPHA = S>*Q 

2 BETA (MIN) =-0.10000 E 01 BET A < M AX ) = 0.10000E 01 

l GAMMA(MIN) =-D.10Q0GE 01 


RU) = 0*0 

PHI = 0.0 


R(Y) = 0.0 
PSI = 0.0 


R(2) = 0.0 

QMEGA= 0.0 REF. TO INTERMEDIATE COORDINATE SYSTEM 15 


INTERMEDIATE COORDINATE SYSTEMS 


ICS NUMBER 15 RX=0.0 


RY=G . 10000E 01 RZ-G * 10 OOOE 01 PHI=0.0 


. F.SlI=0.0^ 



VIEW FACTORS 


element area 


tLEM TO ELEM = VIEW FACTOR 


ELEM TO ELEM - VIEW FACTOR 


VF SUM FROM 
ELEMENT 


100 0 400 0 OE Ol»*»»«»*»***»«»»»*****»»*»****»<>»*««»*****»»*****»»***«**«**?±?*?™?^^^ *t****™ *^^ 

100 - 20 =0*207db0E 00 20 - 100 =0.20T8A9E 00 

iuu cm m + c ** v __ Q »207B5E 00 

n 40000 F n i******************#**«****************«**»*********************i^ ™ 

<iU u « huuuul VI r — 0.20785E 00 


0 

1 

00 



NASA— GSFC 


. J^RQGlRAM SUMMARY __FOR CASE i. 


t PROGRAM TERMINATED NORMALLY 


CPU TIME (SECONDS ) 


___ PROCESSING QE INPUT OAT A= . 0 

VIEW FACTOR COMPUTATI ONS= 0 


NUMBER OF VIEW FACTOR COMPUTATIONS 


FINITE D IFFERENCE= 0 

CONTOUR INTEGRAL= \t 


0 

1 

to 


REQUESTED CORE SPACE ALLOWS THE INPUT DATA GENERATED 


1B0 ELEMENTS 
393 SUB-ELEMENTS 
815 SUB-ELEMENTS 

AFTER COMPRESSION 
5M. GRID, POINTS 


2 ELEMENTS 
8 SUB-ELEMENTS 

4 grid POINTS 


J.6K„cfYTeS..QF. .C0Rt WERE NOT.. USED IN THIS CASE. 


( .36K. ON ELEMENTS# 36K ON SUB-ELEMENTS AND 3bK ON GRID POINTS.) 


<*»#«*#«»# <*»»*« 

?_END__C.ASE_ j 

#**##*#*## 



